D-A040  672 


UNCLASSIFIED 


PENNSYLVANIA  STATE  UNIV  UNIVERSITY  PARK  APPLIED  RESE—ETC  F/6  13/7 
THREE  DIMENSIONAL  INVISCID  EFFECTS  AND  LIMITATIONS  OF  CASCADE  T— ETC(U) 

N00017-73-C-141B 


S 


APR  74  A TAMURA 


TM-74-78 


ML 


4So40672 


AD  A 040  6 


THREE  DIMENSIONAL  INVISCID  EFFECTS  AND  LIMITATIONS 
OF  CASCADE  THEORY  IN  AXIAL  FLOW  TURBOMACHINERY, 


A./Tamura  | 


Technical  Memo rand 
' -PirFe  -No^  TM-74-78 
April  18,  1974 


N#0017-73-C-141 


The  Pennsylvania  State  University 
Institute  for  Science  and  Engineering 
APPLIED  RESEARCH  LABORATORY \/ 

Post  Offj.ce  Box  30 
State  College,  ?A  16801 


NAVY  DEPARTMENT 


NAVAL  ORDNANCE  SYSTEMS  COMMAND 


11 


ACKNOWLEDGMENTS 


The  author  wishes  to  express  sincere  appreciation  to  his  advisor, 
Dr.  B,  Lakshminarayana , for  his  help  and  counsel  during  the  investi- 
gation. The  work  was  sponsored  by  the  Applied  Research  Laboratory  cf 
The  Pennsylvania  State  University  which  operates  under  contract  with 
the  Naval  Ordnance  Systems  Command,  Dr-  M.  T,  Pigott  is  to  be  credited 
for  encouraging  such  a basic  research  which  is  so  vital  in  understanding 
the  complex  flow  phenomena  in  turbomachinery. 


Ill 


TABLE  OF  CONTENTS 


ACKNOWLEDGMENTS  ......................... 

LIST  OF  TABLES  ....... 

LIST  OF  FIGURES  

NOMENCLATURE  .......................... 

I.  INTRODUCTION,  . ............ 

1.1  Statement  of  the  Problem  ............... 

1.2  Review  of  Previous  Investigation  ........... 

1.3  Method  and  Means  of  Investigation,  ........ 

II,  THE  BASIC  THEORY, 

2.1  The  Three-Dimensional  Flow  Field  Due  to  a Single 
Radial  Vortex  Line  of  Constant  Strength  Spanning  an 
Annu lus 

2.2  The  Three-Dimensional  Flow  Field  Due  to  a Single 
Radial  Source  Line  of  Variable  Strength  Spanning  an 

Annulus.  ....................... 

Ill,  APPLICATION  OF  THE  BASIC  THEORY  .............. 

3.1  Superposed  Flow  Field  Due  to  a Series  of  Radial  Vortex 

and  Source  Lines  Replacing  the  Blade  Row  in  Axial 
Turbomachinery  ... 

3-2  Direct  Problem  in  Axial  Turbomachinery  , , , 

IV  NUMERICAL  TECHNIQUE  AND  SOLUTIONS 

4.1  Determination  of  Eigen  Values  and  Related  Functions 
42  Numerical  Calculation  of  the  Three-Dimensional 

Perturbed  Flow  Field  ................ 


4.2.1  Convergency  of  the  Solutions  of  the  Three- 
Dimensional  Perturbed  Flow  Field.  ,...,, 

4.2.2  Three-Dimensional  Perturbed  Flow  Field  Due  to 
a Single  Radial  Vortex  Line  of  Strength  Unity 
and  Its  Comparison  with  Two-Dimensional 
Solution.  .................. 

.3  Three-Dimensional  Perturbed  Flow  Field  Due  to 
a Single  Radial  Source  Line  of  Strength  Unity 
and  Its  Comparison  with  Two-Dimensional 
Solution, 


IV 

TABLE  OF  CONTENTS  (continued) 


Page 


4.3  Numerical  Solution  of  Turbomachinery  Flow  .......  40 

4.3.1  Lifting  Line  Model  ......  40 

4.3.2  Source  Line  Model  43 

4.3.3  Lifting  Surface  Model  Including  Thickness 

Effects.  .....  ....  46 

4.3.4  Effects  of  the  Radial  Source  Lines  of  Variable 

Strength  ........  51 

V,  CONCLUSIONS  .....  .....  54 

REFERENCES 58 

APPENDIX  A:  Specification  of  the  Boundary  Condition  for  Radial 

Vortex  Line  and  Source  Line  59 

APPENDIX  B:  Elimination  of  the  Summation  About  i in  the  Series 

Summation,  , . 62 

APPENDIX  C:  Related  Bessel  Function  and  Characteristics  of  the 

Series  Expression.  64 

APPENDIX  D:  Approximate  Formula  for  the  Solution  of  the  Three- 

Dimensional  Flow  Field  ...............  69 

APPENDIX  E:  Two-Dimensional  Unrolled  Solution,  72 


LIST  OF  TABLES 


Table 

1 Radial  Coordinates  Used  for  Computation  . . . . . 

2 Radial  Coordinates  Used  for  Computation  (v  = 0.6) 

3 Radial  Coordinates  Used  for  Computation  (v  = 0,3) 

4 Strength  of  Vortex  and  Source  Configuration  Used, 


Page 

30 

31 
31 
4? 


VI 


LIST  OF  FIGURES 


Figure 

Page 

1 

Single  Radial  Vortex  and  Source  Line  . 

• • " 

74 

2 

Coordinate  System.  ....... 

. . . 

• • * 

..... 

74 

3 

Blade  Replaced  by  Vortex  and  Source  Line  . , 

..... 

74 

4 

Velocity  Triangle  at  Mid  Radius. 

. . . 

. . , 

4 0 0 0 0 

75 

5a 

H -Constant  Plane  . » 

, . . 

. . . 

75 

5b 

^-Constant  Plane  ........ 

, . . 

. . . 

..... 

75 

6 

Velocity  at  the  Trailing  Edge  to 
Joukowski's  Condition 

Satisfy  Kutta- 

75 

7 

Eigen  Values  for  v = 0,6  , . . . 

. . . 

. . . 

0 0 0 0 0 

76 

8 

Eigen  Values  for  v = 0-3  . . . . 

. . . 

. . . 

..... 

77 

9 

Convergency  of  at  0 - 30°,  r) 

= 0-8 

for  v 

= 0,6  . . 

78 

10 

Convergency  of  usQ  at  Q = 5Qo  ^ p 

= 0.8 

for  v 

= 0.6  . • 

78 

11 

Tangential  Velocity  Distribution 
(0  = 0.6)  Due  to  a Single  Radial 
Strength  Unity  for  v = 0,6  . . . 

uvfl  at 
Vortex 

the  Hub 
Line  of 

79 

12 

Tangential  Velocity  Distribution 
0l  = 0 8)  Due  to  a Single  Radial 

u Qat 
Vortex 

the  Mid  Radius 
Line  of 

Strength  Unity  for  v = 0,6 80 

13  Tangential  Velocity  Distribution  u at  the  Tip 
(n  = 10)  Due  to  a Single  Radial  Vortex  Line  of 

Strength  Unity  for  v = 0,6  ..............  81 

14  Axial  Velocity  Distribution  u at  the  Hub  (o  = 0,6) 

Due  to  a Single  Radial  Vortexv^  Line  of  Strength  Unity 

for  V = 0,6 82 

15  Axial  Velocity  Distribution  u at  the  Mid  Radius 
(n  = 08)  Due  to  a Single  RadYal  Vortex  Line  of 

Strength  Unity  v = 0,6  ........  83 

16  Axial  Velocity  Distribution  u at  the  Tip  (r|  = 1.0) 

Due  to  a Single  Radial  Vortexv  Line  of  Strength  Unity 

for  v*0.6, ........... ..........  84 


Vll 


LIST  OF  FIGURES  (continued) 

Figure  Page 

17  Radial  Velocity  Distribution  u at  the  Mid  Radius 
(r)  = 0-8)  Due  to  a Single  Radial?  Vortex  Line  of 

Strength  Unity  for  V = 0.6,  85 

18  Tangential  Velocity  Distribution  u^at  the  Hub 
(H  ~ 0,6)  Due  to  a Single  Radial  Source  Line  of 

Strength  Unity  for  v = 0,6,  .............  . 86 

19  Tangential  Velocity  Distribution  u q at  the  Mid 
Radius  (n  = 0,8)  Due  to  a Single  Radial  Source  Line 

of  Strength  Unity  for  v=0. 6.. ......  ...  87 

20  Tangential  Velocity  Distribution  u q at  the  Tip 
(H  = 1.0)  Due  to  a Single  Radial  Source  Line  of 

Strength  Unity  for  v = 0,6.  .............  88 

21  Axial  Velocity  Distribution  u at  the  Hub  (n  = 0 6) 

Due  to  a Single  Radial  Sources^Line  of  Strength  Unity 

for  V = 0,6  89 

22  Axial  Velocity  Distribution  u at  the  Mid  Radius 
(n  = 0,8)  Due  to  a Single  Radial  Source  Line  of 

Strength  Unity  for  V*0, 6,  ......  ........  90 

23  Axial  Velocity  Distribution  u at  the  Tip  (n  = 10) 

Due  to  a Single  Radial  SourceS^Line  of  Strength  Unity 

for  v=0. 6.  91 

24  Radial  Velocity  Distribution  u at  the  Mid  Radius 
(H  - 08)  Due  to  a Single  Radill  Source  Line  of 

Strength  Unity  for  V = 0-6  92 

25  Tangential  Velocity  Distribution  u q at  the  Hub 
(n  = 03)  Due  to  a Single  Radial  Vertex  Line  of 

Strength  Unity  for  v =■  0.3  .....  93 

26  Tangential  Velocity  Distribution  u ^ at  the  Mid  Radius 
(r  = 0-65)  Due  to  a Single  Radial  Cortex  Line  of 

Strength  Unity  for  V = 0,3 94 

27  Tangential  Velocity  Distribution  u q at  the  Tip 
(H  = 1.0)  Due  to  a Single  Radial  Vertex  Line  of 

Strength  Unity  for  v = 0 3 95 

28  Axial  Velocity  Distribution  u at  the  Hub  (r|  = 0 3) 

Due  to  a Single  Radial  VortexV^Line  of  Strength  Unity 

for  v » 0 3 . . . , 96 


VI  li 


Figure 

29 


30 


31 


32 


33 


34 


35 


36 


3' 


38 


39 


LIST  OF  FIGURES  (continued) 


Page 

A.xial  Velocity  Distribution  u at  the  Mid  Radius 

(n  = 0,65)  Due  to  a Single  Racial  Vortex  Line  of 

Strength  Unity  foi  V = 0. 3 .......  97 

Axial  Velocity  Distribution  u . at  the  Tip  (n  = 10) 

Due  to  a Single  Radial  Vortexv^Line  of  Strength  Unity 

for  V * 0.3.  . . ...........  98 

Radial  Velocity  Distribution  u „ at  the  Mid  Radius 

(p  = 0,65)  Due  to  a Single  Radiall  Vortex  Line  of 

Strength  Unity  for  V = 0.3  . . . . ..........  99 

Tangential  Velocity  Distribution  u q at  the  Hub 
(P  = 03)  Due  to  a Single  Radial  Source  Line  of 
Strength  Unity  for  V = 0.3  , 100 

Tangential  Velocity  Distribution  u^q  at  the  Mid 

Radius  (p  = 0 65)  Due  to  a Single  Radial  Source  Line 

of  Strength  Unity  for  v = 0,3  101 

Tangential  Velocity  Distribution  u q at  the  Tip 
(p  = 1,0)  Due  to  a Single  Radial  Source  Line  of 
Strength  Unity  for  v = 0,3  ............  - 102 

Axial  Velocity  Distribution  u at  the  Hub  (p  = 03) 

Due  to  a Single  Radial  SourcesHine  of  Strength 

Unity  for  v=  0 3 , . . . , 103 

Axial  Velocity  Distribution  u at  the  Mid  Radius 

(p  - 0 65)  Due  tc  a Single  RaSfal  Source  Line  of 

Strength  Unity  for  v = 0 3 ...........  - . 104 

Axial  Velocity  Distribution  u at  the  Tip  (P  = 1-0) 

Due  to  a Single  Radial  SourceS^Line  of  Strength  Unity 

for  v>  * 0 3 105 

Radial  Velocity  Distribution  u n at  the  Mid  Radius 

(n  = 0 65)  Due  to  a Single  Radfal  Source  Line  of 

Strength  Unity  for  v * 0.3  , 106 

Ratio  of  the  Tangential  Velocities  Induced  by  the 

Radial  Vortex  Lines  of  JVS  = 1 at  0 = 2°  107 


40 


Ratio  of  the  Tangential  Velocities  Induced  by  the 
Radial  Vortex  Lines  of  JVS  = 1 at  0 = 9°  . , . , 


108 


IX 


LIST  OF  FIGURES  (continued) 

Figure 

Page 

41 

Ratio  of  the  Tangential  Velocities  Induced  by  the 
Radial  Vortex  Lines  of  JVS  = 1 at  0 = 1/2,  1/4,  1/8, 
1/16-Spacing  Points  for  NB  = 3 ...........  . 

109 

42a  ,b 

Ratio  of  the  Axial  Velocities  Induced  by  the  Radial 
Vortex  Lines  of  JVS  = 1 at  £ = 0.1,  n = v and  l + v/2 

110 

42c 

Ratio  of  the  Axial  Velocities  Induced  by  the  Radial 
Vortex  Lines  of  JVS  =1  at  £ = 0,1,  r|  = 1 0 - , . , . 

111 

43 

Ratio  of  the  Radial  Velocities  Induced  by  the  Radial 
Vortex  Lines  of  JVS  = 1 at  0 = 2°,  p = l+v/2  . . . , , 

112 

44 

Ratio  of  the  Radial  Velocities  Induced  by  the  Radial 
Vortex  Lines  of  JVS  = 1 at  0 = 9°,  r|  = l+v/2  . . « . , 

112 

45 

Ratio  of  the  Tangential  Velocities  Induced  by  the 
Radial  Source  Lines  of  JVS  =lat£=0.1,n=v, 

l+V/2  and  1.0.  ...................  . 

1 13 

46 

Ratio  of  the  Axial  Velocities  Induced  by  the  Radial 
Source  Lines  of  JVS  = 1 at  0 = 2°,  0 = v , l+v/2  and 

1,0.  ......  

114 

47 

Ratic  of  the  Radial  Velocities  Induced  by  the  Radial 
Source  Lines  cf  JVS  = 1 at  0 = 2°,  p = l+v/2  

115 

48 

Sketch  of  a Model  Including  Eight  Radial  Vortex  and 
Source  Lines 

116 

49 

Singular  Points  and  Calculation  Points  on  (£',  n . ,0)- 
m raid 

117 

50 

Singular  Points  and  Calculation  Points  on  (£',  0) - 
Plane 

118 

51 

Distribution  of  the  Tangential  Velocities  at  the  Hub 
Induced  by  NB  Series  of  Radial  Vortex  Lines  of  JVS  = 8 
for  A = 0 , NB  = 1 and  V=0. 3.. .......... 

119 

52 

Distribution  of  the  Tangential  Velocities  at  the  Hub 

Induced  by  NB  Series  of  Radial  Vortex  Lines  of  JVS  = 8 

for  A = 54°,  NB  = 1 and  v = 0.3  120 


53  Ratio  of  the  Tangential  Velocities  Induced  by  NB 

Series  of  Radial  Vortex  Lines  of  JVS  * 8 for  A9  = 2° 
at  n = v , l+v/2  and  1.0  , 


121 


Figure. 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 


LIST  OF  FIGURES  (continued) 


Distribution  of  the  Tangential  Velocities  at  the  Hub 
Induced  by  NB  Series  of  Radial  Vortex  Lines  of  JVS  = 

8 for  A = 18°,  NB  = 3 and  V = 0.3  ■ 

Distribution  of  the  Tangential  Velocities  at  the  Tip 
Induced  bv  NB  Series  of  Radial  Vortex  Lines  of  JVS  = 

8 for  A =’ 18  , NB  * 3 and  v = 0.3  

Distribution  of  the  Axial  Velocities  at  the  Hub  Induced 
by  NB  Series  of  Radial  Vortex  Lines  of  JVS  = 8 for  A = 

0 , NB  = 1 and  v = 0.3 . . . . . . . . . , 

Distribution  of  the  Axial  Velocities  at  the  Hub  Induced 
by  NB  Series  of  Radial  Vortex  Lines  of  JVS  = 8 for 
A = 54°,  NB  = 1 and  * 0.3  .............  . 

Ratio  of  the  Axial  Velocities  at  the  Hub  Induced  by  NB 
Series  cf  Radial  Vortex  Lines  of  JVS  - 8 for  NB  = 1 at 
C/c  cos  A=  0 0 

Distribution  of  the  Radial  Velocities  at  the  Mid  Radius 
Induced  by  NB  Series  of  Radial  Vortex  Lines  of  JVS  = 8 

for  A = 18°  at  C/c  cos  A=  0.0  ............. 

Ratio  of  the  Tangential  Velocities  at  the  Hub  Induced 
by  NB  Series  of  Radial  Source  Lines  of  JVS  = 8 for  NB  = 

1 at  C/c  cos  A = 0,0  . . ........  ....... 

Ratio  of  the  Tangential  Vel  'ities  at  the  Hub  Induced 

by  NB  Series  of  Radial  Source  Lines  of  JVS  = 8 for  NB  = 

1 at  ^/c  cos  A = 0.0  , . . . , . . . . . . . . « , . , 

Ratio  of  the  Axial  Velocities  Induced  by  NB  Series  of 

Radial  Source  Lines  of  JVS  = 8 for  A = 0°  at  A0  = 2 , 

1 = V,  l+v/2,  and  1.0  ................  , 

Ratio  of  the  Axial  Velocities  at  the  Hub  Induced  by  NB 
Series  of  Radial  Source  Lines  of  JVS  = 8 for  NB  = 20  at 


Distribution  of  the  Radial  Velocities  at  the  Mid  Radius 
Induced  by  NB  Series  of  Radial  Source  Lines  of  JVS  = 8 
for  A = 18°  at  £/c  cos  A = 0,0  ........... 

Distribution  of  the  Source  of  Variable  Strength  and 
Coordinate  of  the  Singular  Points  and  Calculation 
Points 


x 


Page 


122 


123 


124 


125 


126 


1 27 


128 


] 28 


129 


130 


130 


131 


XI 


LIST  OF  FIGURES  (continued) 


Figure  Page 


66 

Pressure  Distribution  Due  to  Source 
Variable  Strength  at  £ = 0=0  for  NB 

Lines  of 

= 3 and  V = 0.6  - 

132 

67 

Pressure  Distribution  Due  to  Source 
Variable  Strength  at  £ = 0=5  for  NB 

Lines  of 

= 3 and  v = 0 6 . 

132 

68 

Pressure  Distribution  Due  to  Source 
Constant  Strength  Qj.  i = 0,2,  Qi  2 = 
£ = 0.0  for  NB  = 3 and  v = 0.6.  ! » 

Lines  of 
- 0 2 at 

133 

69 

Pressure  Distribution  Due  to  Source 
Constant  Strength  = 0.2,  Qi  2 = 

X,  - 0,5  for  NB  = 3 and  v = 0.6  ’ , 

Lines  of 
- 0,2  at 

133 

70 

Split  Radial  Vortex  Line.  ..... 

1 0 4 

71 

Split  Radial  Source  Line.  ..... 

* , t e <*  » . 

. 134 

Xil 


NOMENCLATURE 


Rh 

0 


Hub-to-tip  ratio 

Radius  of  the  outer  casing  wall 

Radius  of  the  hub 

Angular  coordinate  in  a cylindrical  coordinate  system 
(rad) 

Tangential  coordinate  of  the  j -th  singular  point  on  i-rh 
blade 


n Nondimensional  radial  coordinate  normalized  by  R^  in  a 

cylindrical  coordinate  system 

H . , Mid  radius  = (1  + v)/2 

mid 

C Nondimensional  axial  coordinate  normalized  by  R in  a 

cylindrical  coordinate  system 


C . 
'3 


Nondimensional  axial  coordinate  of  the  j -th  singular 
point 


r. 

3 


Qj  (") 


0 


Nondimensional  strength  of  the  j -th  radial  vortex  line 

normalized  by  R.  • W_ 

’ t Cm 

Nondimensional  strength  of  the  j -th  radial  source  line 

normalized  by  R_  • W, 
t Cm 

Strength  of  the  split  vortex  line 


p Strength  of  the  split  source  line 

6 Split  vortex  and  source  line  distributed  over  j 0 j £ 6/2 


$ 

v 


$ 

s 


uv0,uvc’uvn 


Nondimensional  velocity  potential  due  to  a single  radial 
vortex  line  or  series  of  radial  vortex  lines  of  constant 
strength  normalized  by  R • IV 

Nondimensional  velocity  potential  due  to  a single  radial 

source  line  or  series  of  radial  source  lines  of  variable 

strength  normalized  by  R^  • W, 

6 1 t Cm 

Nondimensional  tangential,  radial  and  axial  velocity 
components  due  to  a single  radial  vortex  line  or  series 
of  radial  vortex  lines  of  constant  strength 


| ^ | M | 


J 


NOMENCLATURE  (continued) 


Nondimensional  tangential,  radial  and  axial  velocity 
components  due  to  a single  radial  source  line  or  series 
of  radial  source  lines  of  variable  strength 


i,  j,  k 


Unit  vector  in  0,  r)  and  £ directions 


WWW 

0’  V \ 


(a  ) 
mm 


(a. ) 

1 m 


V V V 

0’  V t 


Induced  velocity  vector  due  to  a series  of  radial  vortex 
lines  of  constant  strength 

Induced  velocity  vector  due  to  a series  of  radial  source 
lines  of  variable  strength 

Three-dimensional  cascade  mean  velocity  vector 
Tangential,  radial  and  axial  components  of  W 

(Vn  = l+v/2 

(lVn  = l+v/2 
Cascade  mean  flow  angle 

iam'Jr|  = l+v/2 
Cascade  inlet  flow  angle 

^°i^n  = l+v/2 
Resultant  velocity  vector 

Tangential,  radial  and  axial  components  of  V 
Number  of  singular  points  per  blade 
Number  of  blades 

Order  of  the  polynomial  for  strength  of  the  source  Q(j.p) 
Stagger  angle 


= (A) 


n = l+v/2 


Spacing  of  the  blade  row 
Chord  length  of  the  blade 


Solidity  of  the  blade  row 


Pressure  coefficient 


XIV 


mn 

J (K  n) 
m mn 

Y (K  n) 
m'  mn  ' 

Z (K  ,n) 
m mn 


F (K  ,v) 
v mn 


F (K  ,v,j) 
s mn 


SUBSCRIPTS 


3 

v 

s 

C 


NOMENCLATURE  (continued) 

n-th  eigen  value  referred  to  m-th  order  Bessel  function 

The  first  kind  Bessel  function  of  m-th  order 

The  second  kind  Bessel  function  of  m-th  order 

= J (K  n)  * Y ' (K  ) - J ' (K  ) • Y (K  n) 
m mn  m mn  m mn  m mn  ' 


= f1  nz  (k  n)  dn/  f1  n zj  (k  n)  dn 

4 m mn  4 ™ 


v 

1 


m mn 
1 


= / o.(n)z  (k  n)  dn/  / n z (k  n)  dn 

v m mn  J ^ m ^ mn  J 


Refer  to  m-th  order  Bessel  function 

Refer  to  n-th  eigen  value  for  each  order  of  the  Bessel 
function 

Refer  to  i-th  blade 

Refer  to  j-th  singular  point  for  any  blade 
Refer  to  the  vortex  system 
Refer  to  the  source  system 
Refer  to  mid  radius 


CHAPTER  I 


INTRODUCTION 

1 . 1 Statement  of  the  Problem 

Axial -flow-propeller  pumps,  inducers,  or  blowers  consist  of  a few 
blades  of  high  solidity  (chord  length-to-spacing  ratio)  The  chord 
length  of  these  blades  is  comparable  to  the  hub  diameter,  and  the  blades 
usually  overlap,  As  a result  of  these  characteristics  and  the  low  hub 
to-tip  ratio  usually  encountered  in  practice,  the  flow  field  through 
such  an  impeller  is  three-dimensional  In  order  to  solve  this  kind  of 
flow  field,  two-dimensional  unrolled  potential  flow  theory  (which  is 
two-dimensional  potential  flow  theory  adapted  to  a cylindrical 
coordinate  system)  or  streamline  curvature  methods  are  available.  The 
two-dimensional  unrolled  solution  gives  tangential  and  axial  velocity 
components  only  Furthermore,  there  is  no  interference  effect  between 
each  radius,  especially  hub  and  tip,  since  the  two-dimensional  unrolled 
solution  is  applied  to  each  radius  independently.  Alternatively,  the 
streamline  curvature  method,  which  is  the  most  powerful  method  in  the 
analysis  of  turbomachinery,  is  essentially  the  theory  for  axisymmetric 
flow  Thus,  it  is  not  suitable  for  evaluating  the  local  three- 
dimensional  flow  field  because  of  the  few  number  of  blades  The  stream 
line  curvature  method  is  very  powerful  for  obtaining  the  correlation  of 
inlet  and  exit  flow  when  the  flow  field  is  axisymmetric.  To  solve  the 
three-dimensional  direct  problem  in  turbomachinery  (that  is,  to  deter- 
mine the  flow  field  around  a given  blade  profile)  by  the  relaxation 
method  will  be  extremely  time  consuming  since  the  boundary  shape  is 


r 
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complicated.  Therefore,  no  proper  method  exists  for  the  design  or 
analysis  of  the  three-dimensional  flow  through  such  an  impeller. 


1 . 2 Review  of  Previous  Investigation 

The  basic  analysis  of  the  three-dimensional  potential  flow  field 
in  axial  flow  turbomachinery  carried  out  by  Meyer  (1)  was  an  extension 
of  Ackeret's  idea  for  two-dimensional  potential  flow  around  a cascade. 
The  theoretical  model  was  made  up  of  ring  and  line  vortices  to  replace 
the  blades  and  boundary  wall, 

Tyson  (2)  showed  that  the  discrepancy  between  the  three-dimensional 
and  the  two-dimensional  unrolled  solutions  is  large  for  a radial  vortex 
line  spanning  an  annulus,  A single  radial  vortex  line  of  constant 
strength  spans  an  annulus  bounded  by  a hub  and  outer  casing  wall-  The 
whole  flow  field  except  the  radial  vortex  line  is  irrotational  The 
basic  equation  reduces  to  Laplace's  equation  and  solved  with  the 
following  boundary  conditions  (see  Figures  1 and  2): 

(a)  Radial  velocity  at  hub  and  outer  casing  wall  must  vanish. 

(b)  Axial  and  radial  velocity  must  vanish  at  £ = ± 

(c)  Tangential  velocity  at  £ = ± 00  must  be  ±r/4nri  to  satisfy 
irrotationality . 

(d)  Boundary  conditions  were  specified  by 


1_ 

n 
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£->±o 


The  numerical  results  obtained  by  Tyson  showed  100%  difference  between 
the  three-dimensional  and  two-dimensional  unrolled  solutions  in  both 
tangential  and  axial  velocity  components  at  hub  and  outer  casing  wal1 
for  hub -to -tip  ratio  v = 0,6,  but  they  coincided  with  the  two- 
dimensional  solution  at  mid-radius.  However,  Etter,  Van  Dyke  and  the 
present  author  obtained  different  results.  The  same  calculations  showed 
that  the  discrepancy  between  two-dimensional  unrolled  solution  and  the 
three-dimensional  solution  is  less  than  those  calculated  by  Tyson. 

Rossow  (3)  developed  the  three-dimensional  potential  theory  for  a 
single  radial  vortex  line  of  constant  strength  and  a single  radial 
source  line  of  constant  strength  spanning  an  annulus.  The  analysis  for 
the  radial  vortex  line  was  essentially  the  same  as  that  of  Tyson's,  but 
the  results  were  given  in  terms  of  Lommel's  function,  A solution  for  a 
radial  source  line  of  constant  strength  was  developed  similar  to  that  of 
the  radial  vortex. 

Etter  and  Van  Dyke  (4)  calculated  the  three-dimensional  flow  field 
due  to  a radial  vortex  line  in  a cylindrical  annulus  based  on  Tyson's 
analysis.  Eigen  values  and  related  functions  for  hub-to-tip  ratios 
v = 0,6,  07  and  08  were  obtained  for  Bessel  functions  of  orders  m = 1 
to  40,  the  number  of  eigen  values  foi  each  fixed  value  of  m is  n = 1 to 
6 The  velocity  distribution  due  to  a single  radial  vortex  line  of 
constant  strength  was  different  from  Tyson's  (2),  but  coincident  with 
the  author's  results.  Example  calculations  were  made  for  a four-bladed 
cascade  using  6 vortex  lines  per  blade  for  hub-to-tip  ratio  v = 0.7, 

The  discrepancies  between  the  two-dimensional  solution  and  three- 
dimensional  solution  were  quite  small  for  the  mean  tangential  velocity 


distribution 
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1 . 3 Method  and  Means  of  Investigation 

The  problem  is  to  develop  a theory  for  predicting  the  three- 
dimensional  potential  flow  in  axial  flow  turbomachinery.  The  general 
objective  of  the  investigation  will  be  to  obtain  a reliable  understanding 
of  the  three-dimensional  potential  flow  through  axial  flow  turbomachinery 
The  calculation  of  the  three-dimensional  potential  flow  through  an 
impeller  is  based  on  the  method  of  distributed  singularities.  The 
blades  are  replaced  by  a series  of  line  vortices  and  line  sources  which 
have  their  axes  along  the  radial  direction  and  are  arranged  along  the 
blade  camber  surface.  The  basic  perturbed  velocity  fields  due  to  a 
single  radial  vortex  line  of  constant  strength  and  a single  radial 
source  line  of  variable  strength  along  the  radial  direction  were 
computed  from  a modified  theory  based  on  Tyson's  and  Rossow's  formula- 
tions. Since  the  resultant  perturbed  flow  fields  are  solutions  of  the 
linear  Laplace  equation  and  the  boundary  conditions  are  homogeneous, 
these  solutions  may  be  superimposed  for  a series  of  radial  vortex  lines 
and  source  lines  located  at  arbitrary  points.  The  perturbed  flow  field 
due  to  a single  radial  vortex  line  of  constant  strength  was  compared 
with  Tyson's  and  Etter's  results  for  hub-to-tip  ratio  v = 06.  The 
results  are  coincident  with  Etter's  but  different  from  Tyson's  The 
perturbed  flow  field  due  to  a single  radial  source  line  of  constant 
strength  shows  good  agreement  with  the  two-dimensional  solution  at  the 
mid  radius.  Eigen  values  and  related  functions  were  calculated  for 
Bessel  functions  of  orders  from  0 to  100,  the  number  of  eigen  values  for 
each  order  of  Bessel  function  n = 20,  hub-to-tip  ratio  v = 03,  06 
The  number  of  eigen  values  and  the  order  of  the  Bessel  functions  required 
to  obtain  enough  accuracy  increase  greatly  in  the  vicinity  of  the 


singular  points.  Also,  for  smaller  hub-to-tip  ratios,  more  eigen  values 
are  required.  A computer  program  to  obtain  the  three-dimensional  flow 
field  for  an  arbitrary  number  of  blades,  number  of  singular  points  per 
blade,  and  singular  point  locations  was  developed.  Examples  exhibiting 
the  interference  effects  due  to  hub-to-tip  ratio,  stagger  angle,  and 
number  of  the  blades  were  carried  out.  The  effects  of  the  radial 
variation  of  the  strength  of  the  radial  source  line  were  also  examined 
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CHAPTER  II 
THE  BASIC  THEORY 


2.1  The  Three-Dimensional  Flow  Field  Due  to 


a Single  Radial  Vortex  Line  of  Constant 
Strength  Spanning  an  Annulus 


In  order  to  carry  out  a potential  flow  analysis,  the  following 
basic  assumptions  were  made:  the  flow  is  inviscid,  incompressible  and 

has  no  heat  transfer.  A single  radial  vortex  line  of  constant  strength 
F spans  the  annulus  region  between  the  hub  and  outer  casing  wall  (see 
Figure  1) . The  theoretical  work  reported  in  this  section  is  due  to 
Tyson  (2) . The  whole  flow  field  is  irrotational  except  for  the  vortex 
line.  The  irrotational  and  continuity  conditions  are  reduced  to  a 


Laplace  equation  in  terms  of  the  velocity  potential  and  0,  n,  C,  (see 


Figure  1) , 


923> 


9 1 2 1 

n2  . _ 1 9 ^ V ..  1 V v 

(n3n 


= 0 


(1) 


The  boundary  conditions  are: 


a)  The  axial  velocity  should  vanish  at  £ 

9$ 


(ar} 


= 0 


£-*±a 


(2) 


b)  The  radial  velocity  should  vanish  at  hub  and  outer  casing  wall 

9$ 

= 0 (3) 

n=v,l 


c)  The  radial  velocity  should  vanish  at  £ = ±°°, 

9$ 


<5T> 


= 0 


(4) 
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d)  At  plus  and  minus  infinity,  the  distribution  of  the  tangential 
velocity  must  be  irrotational , Furthermore,  the  total  magnitude  of  the 
circulation  must  be  equal  to  the  strength  of  the  radial  vortex  line. 

For  convenience,  let  the  magnitude  of  the  circulation  around  the  hub  on 
either  side  of  the  vortex  line  be  17 2,  These  boundary  conditions  are 
the  same  as  those  of  Tyson  (Reference  2)  but  the  definition  of  the 
positive  direction  of  rotation  of  the  vortex  is  based  on  the  right-hand 
screw  rule.  The  corresponding  boundary  condition  of  Reference  (3) 
appears  different  from  the  above  expression,  but  both  are  quite  true 
The  only  difference  is  that  the  initial  value  of  the  circulation 
upstream  due  to  the  vortex  line  is  taken  to  be  -T/2  here  and  zero  in 
Reference  3.  The  present  objective  is  to  determine  the  perturbed  flow 
field.  The  above  difference  of  initial  circulation  may  easily  be  taken 
into  consideration  by  adjusting  the  mean  flow  condition  which  will  be 
superposed  on  the  above  perturbed  flow.  Circulation  about  the  £ axis 
is  -T/2  at  £ = + 00  and  T/2  at  £ = - then,  the  tangential  velocity  at 
£ = ±oo  is  given  by  the  following  equations: 


, 34> 

n l90  J 


£->■+00 


r 

4nn 


(5) 


and 

9$  r 

n(5e%^,=  * (6) 

e)  The  boundary  conditions  must  be  specified  in  the  vicinity  of 
the  vortex  line.  The  behavior  of  the  tangential  velocity  close  to  the 
singular  point  (£  = O-plane)  is  represented  by  the  following  expression; 
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. 3$ 

n '•ae  \ 


}<*<7 


|<  e<2^4 


Detailed  derivation  of  the  specification  of  the  boundary  conditions  is 
given  in  Appendix  A,  The  solution  of  the  Laplace  equation  is  given  by 
the  following  expression: 


$ = I (A  cos  m9  + B sin  m0)  E (C*e  + D*e^)  (aJ  (Kr|)  + b Y (Kn)), 

v . m m ' , nr  J m J J ’ 

m=o  n=l 


where  a,b,  A , B^,  C,  D,  K are  constants  which  must  be  determined,  m,  n 

are  integer  numbers  and  J , Y are  the  first  and  second  kind  of  m-th 

mm 

order  Bessel  functions.  From  boundary  condition  (a),  D = 0 when  £ > 0, 
and  C = 0 when  £ < 0.  Without  loss  of  generality,  C and  D can  be  taken 
as  unity.  Using  boundary  condition  (b) , substitution  of  Equation  (8) 
into  Equation  (3)  gives  the  following  conditions: 


a J ' (Kv)  + b Y ' (Kv)  = 0 
m ' m ' 


where 


a J ' (K)  + b Y ' (K)  =0 

m v ' m v ' 


3J  (Kn) 


V ™ ■ -3?r 


3Ym(Kn) 


* V ™ - “5T 


Combining  Equations  (9)  and  (10), 


J ' 
m 

(K  ) 
mn' 

J ' 
m 

(K  v) 
mn 

Y ' 
m 

(K  ) 
mn' 

Y ' 
m 

(K  V) 
v mn  ’ 

(ID 
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Equation  (11)  is  a relationship  between  the  hub-to-tip  ratio  v and  the 
eigen  value  Kmn,  where  the  subscript  m refers  to  the  m-th  order  Bessel 
function  and  n is  the  n-th  eigen  value  for  each  m,  (m  = 1,2,3,,..,  n = 1, 
2,3,....  ).  Introducing  a new  expression  zm(Kmnn) > the  last  term  of 
Equation  (8)  is  given  by  the  equation. 


a J (K  n)  + b Y (K  n)  = 2 Z (K  n)  , (12) 

nr  mn  ' nr  mn  ' Y ' (K  ) m mn  ' * v 1 

m ^ nur 

where 


Z (K  n)  = J (K  n)  * Y ' (K  ) - J ' (K  )•  Y (K  n)  • 
nr  mn  ' nr  mn  ' m ^ nur  m v mn-1  m mn 


Therefore,  Equation  (8)  is  reduced  to  the  following  expression: 


4>  = E (A  cos  m 0 + B sin  m0)  E exp  (+K  £,)  Z (K  q) 

V m=o  m m n=l  mn  Y ' (K  p)  m 1,111 

m mn 


The  relation  A , B , and  a are  determined  by  the  orthogonal  charac- 
m m ° 

teristics  of  the  trigonometric  functions  and  Bessel  function.  In  order 
to  apply  the  boundary  condition  (e),  differentiate  with  respect  to 
0 and  equate  it  with  Equation  (7) , 

00  oo 

[I  E m (-A  sin  m0  ♦ Bm  (cos  m0)  E ,a  - y W)  ] 

1 m=o  n=l  m v mn  £-*±o 


-I  <»<  I 


7 < 9 < - r 
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Multiply  both  sides  by  (cos  m 0 + i sin  m 0)  • n (K  n) , then 
integrate  from  - tt  to  tt,  and  from  v to  1 for  p: 


m<-\  - 1A»>  yfiry  I dr>  ■ ; 4/'  dn 


m 2 'v  " rnn' 


Taking  the  limit  6 ■+■  0,  and  holding  6y  = T constant,  then. 


(Bm  * lAm}  Y ' (K  ) ^ nZm  (Kmnn)  dn  = + 2 ^ nZm(Kmnr|)  dn  ‘ 


m mn'  v 


From  the  imaginary  part, 


From  the  real  part. 


A = 0 . 
m 


B = 
m 2imt 


/ nz  (K  n)dn 

' v m mn 


Y ' (K  ) 71  ~2rv  , 

m m*  / nzm  (Kmnn)  dn 


where  m=l,2,3,.  . . . 

Taking  into  consideration  the  boundary  condition  (b) , the  solution  of 
Laplace  Equation  (1)  is  given  by  Equation  (20): 


; (0  + 2 E,  5±IL“L  z exp  (-K  c)  F (K  v)  Z (K  n)), 

4tt  m=l  m n_^  r v+  mnw  mn,  m mn  ; 1 ' 


J n Z (K  n)  dn 

' m mn 

F (K  ,v)  = -r = 

v 11111  r i 2 rv 

I n z (K  n)  dn 

; m ^ mnu 


Therefore,  the  tangential,  radial,  and  axial  velocity  components  are 
given  by  the  following  expressions: 
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u . = ; -jd—  (1  + 2 I cos  m0  E exp  (ZK  £)  • F (K  ,v)  • Z (K  n)), 
v0  4irn  . 1 mn^'  vv  mn’  m mn 

m=l  n=l 

(22) 

r 00  . 0°  3Z  (K  n) 

— l n sin  mo  » f —1/  _ r.  fv  % m mn  , » 

u = + =—  E E exp  (+K  £)  • F (K  ,v)  • r (23) 

vn  2tt  . m •,  mn  v mn’  ^ 3n  v J 

m=l  n=l 

and 


u = ± - E S1'n~m9'  E K exp  (;K  5)  • F (K  ,v)  • Z (K  n)  , 

vC  2tt  , m , mn  r mn^J  v v mn’  mv  mn  J * 

m=1  n=1  (24) 


where  the  upper  and  lower  signs  refer  to  +£  and  -£  directions 
respectively. 


2.2  The  Three-Dimensional  Flow  Field  Due 

to  a Single  Radial  Source  Line  of  Variable 
Strength  Spanning  an  Annulus 

The  basic  assumptions  regarding  the  flow  field  (irrotational , 
inviscid,  no  heat  transfer,  and  incompressible)  are  the  same  as  those 
listed  in  che  previous  section.  This  section  is  based  on  Rossow's 
analysis  (Reference  3),  The  mathematical  model  of  a single  radial 
source  line  is  similar  to  that  of  a single  radial  vortex  line  except 
that  the  radial  source  line  may  have  a variable  strength  along  the 
radial  direction.  In  the  case  of  a single  radial  vortex  line,  if  the 
strength  of  the  vortex  varies  in  the  radial  direction,  there  will  be 
vortices  shed  according  to  Helmholtz's  law.  The  flow  field  is  no 
longer  irrotational,  but  the  flow  field  due  to  a variable  strength 
radial  source  line  is  still  irrotational.  A single  radial  source  line 
of  variable  strength  Q(n)  spans  the  annulus  region  between  the  hub  and 

L . A 
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outer  casing  wall  (see  Figure  1).  As  in  the  previous  section,  the 
whole  flow  field  is  irrotational , and  the  velocity  potential  exists, 
satisfying  the  Laplace  equation.  Let  the  velocity  potential  due  to  a 


single  radial  source  line  be  . Therefore, 

. * 3$  . 32$  32<i> 

FT  w ">  7 ^2-  * ^2-  ■ • < ) 

The  boundary  conditions  are: 

a)  One -half  of  the  mass  flow  due  to  a single  radial  source  line 
goes  to  £ = +00  and  the  other  half  goes  to  ( = Therefore,  from  the 

continuity  relationship. 


± / Q(n)  dn/2  = (jj-1)  • tt  (l  -v2)  . 


b)  The  radial  velocity  should  vanish  at  hub  and  outer  casing  wall, 


far*  = 0 


n=i,v 


c)  The  radial  velocity  should  vanish  at  £ = ±°°, 


(3Tr)  = 0 


d)  The  tangential  velocity  should  vanish  at  ( = ±°°, 


. 3$ 

- ° ' 


13 


e)  The  specification  of  the  condition  at  points  close  to  the 
single  radial  source  line  (£  = 0-plane)  is  given  by  Equation  (30)  The 


detailed  derivation  is  given  in  Appendix  A. 

( ° • I < 6 < 2"  - r 


<sr> 


; u(n)  6 < a < 6 

Hn”  ? 0 2 


where  the  upper  and  lower  signs  refer  to  £ > 0 and  £ < 0 respectively. 
The  solution  of  the  Laplace  equation  under  the  given  boundary  conditions 
is  given  by  the  following  equation,  which  is  similar  to  those  given  m 
the  previous  section: 


1 00 

*s  = 57  7~T  f ^ dn  - n=2  exP  (+_Kon?) 


F (K  ,v) 
s on 


1 -v  v 


Z (K  r) 
o on 


where 


oo  oo  F (K  ,V) 

-2  Z cos  m6  Z exp  (+K  £)  • — m- • Z (K  n) ) , 

, . r K mv  mn 

m=  1 n=  1 mn 


/„  dn 


F (K  ,v)  = -i t 

s mn’  fl  2 , , 

n z (K  n)  dn 

J m mn 

v 


The  tangential,  radial  and  axial  velocity  components  due  to  a single 
radial  source  line  are  given  by  the  following  expressions, 


°°  00  F (K  ,v) 

u - = £ m sin  m0  Z exp  (+K  C)  -tt-122 * ZfK  n)  , (33) 

s6  2-rrn  m=1  n=1  mn  Kmn  m mn 
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CHAPTER  III 

APPLICATION  OF  THE  BASIC  THEORY 


3 , 1 Superposed  Flow  Field  Due  to  a Series 

of  Radial  Vortex  and  Source  Lines  Replacing 
the  Blade  Row  in  Axial  Turbomachinery 

In  the  preceding  paragraph,  the  three-dimensional  perturbed  flow 

field  due  to  a single  radial  vortex  and  source  line  were  discussed- 

These  basic  solutions  will  be  used  to  construct  a three-dimensional 

theory  to  evaluate  the  potential  flow  field  in  a turbomachinery  blade 

row.  In  general,  the  three-dimensional  flow  field  around  a blade  row 

can  be  calculated  by  means  of  superposition  of  the  three -dimens. onal 

cascade  mean  flow  and  perturbed  flow  field.  The  three-dimensional 

perturbed  flow  field  is  constructed  by  the  superposition  of  the  flow 

fields  due  to  a series  of  radial  vortex  lines  of  constant  strength  and 

source  lines  of  variable  strength  by  which  the  blades  are  replaced- 

This  is  valid  for  a free  vortex  type  blade  row  with  any  arbitrary  blade 

thickness  The  three-dimensional  cascade  mean  flow  is  similar  to  the 

two-dimensional  flow,  but  varying  along  the  radia1  direction  Also, 

strictly  speaking,  the  distribution  of  the  cascade  mean  velocities  must 

be  irrotational , Let  W be  a three-dimensional  cascade  mean  velocity 

observed  from  the  coordinate  system  fixed  to  the  blade  and  WQ,  W and 

o rj 

are  the  components  in  the  0,  n,  £ directions  respectively.  Then, 
for  a free  vortex  blade  row. 


p • W = constant, 

U 


W =0 
n 


(36) 

(37) 


and 


constant  , 


(38) 
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However,  if  the  three-dimensional  perturbed  flow  field  is  small  compared 
to  the  three-dimensional  cascade  mean  velocity,  it  seems  that  the 
restriction  on  the  cascade  mean  velocity  given  by  Equation  (36)  to  (38) 
is  not  compulsarily  needed.  All  the  solutions  for  the  three- 
dimensional  perturbed  flow  field  are  solutions  of  the  linear  Laplace 
equation  and  the  boundary  conditions  are  homogeneous  Therefore, 
mathematically,  the  method  of  superposition  is  valid.  Hence,  the 
resultant  three-dimensional  flow  field  around  a blade  row  will  be 
determined  by  the  superposition  of  the  three-dimensional  cascade  mean 
flow  and  the  perturbed  flow  due  to  both  radial  vortex  lines  and  source 
lines  of  known  magnitude.  Let  4>v  be  the  superposed  three-dimensional 
perturbed  velocity  potential  due  to  a series  of  radial  vortex  lines 
(of  constant  strength)  replacing  a row  of  equally  spaced  blades  (Figure  3) 


JVS 

$ = + £ 

V J-l 


r . NB  00  sin  m(0  - 0 . . ) 

£ (0  - 0.  . + 2 £ — 

4tt  . , 13  , 

i=l  m=l 


m 


£ exp  (4  K (5  - ?.))  • F (K  ,v)  • Z (K  n) ) , (39) 

mn  * j v mn  m mn  J ’ 1 


n=l 


where  NB  is  the  number  of  blades,  0^  is  the  angular  coordinate  for 

j -th  radial  vortex  line  of  the  i-th  blade,  £.  is  the  axial  coordinate 
j J 

of  the  j -th  radial  vortex  line  for  any  blade,  and  I\  is  the  strength  of 
the  3 -th  vortex  line  for  any  blade.  Since  the  blades  are  equally  spaced 
and  have  the  profile,  the  distribution  of  the  vortex  lines  is  the  same 
for  each  blade.  Then,  the  .angular  coordinate  corresponding  to  3-th 
vortex  line  of  i-th  blade  0^,^  is  periodic.  Hence, 


0.  . = 0 . ♦ (i  - 1)  . 

13  I3  NB 


(40) 
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Using  Equation  (40),  the  summation  with  respect  to  i will  be  eliminated 
mathematically,  and  also  the  summation  with  respect  to  m remains  only 
when  m is  a multiple  of  the  number  of  the  blade  NB.  Detailed  derivation 
is  given  in  Appendix  B.  Therefore,  Equation  (39)  is  rewritten  as 
follows,  neglecting  constant  terms.  For  simplicity,  let 


Then, 


k = l • NB 


(41) 


4> 

v 


4"  j.i 


r.  (6  + 2 e 

3 JUi 


sin  k (0  - 01  .) 

•I  9 3 

k 


' = ex?  <♦-  Kkn  « - F/Kkn’V>  ^ 

n=l  J 

A closed  form  of  the  series  expression  for  the  other  summation  has  not 
been  found  so  far,  but  characteristics  of  the  series  expression  with 
respect  to  the  Bessel  function  is  given  in  Appendix  C,  The  corresponding 
three-dimensional  perturbed  flow  field  due  to  a series  of  radial  vortex 
lines  of  constant  strength  replacing  an  equally  spaced  blade  row  of  NB - 
blades  is  given  by  the  following  expressions. 


u 


v0 


. NB_ 
+ 4irq 


JVS 

E r . (1+2  E cos  k (0  - 0 ) 

j=l  3 £=1  3 


Z exp  («  « - 5 ))  • F (Kkn,v)  Z (K  n))  . 

n=l  J 


(43) 


potential  due  to  a series  of  radial  source  lines  of  variable  strength 
replacing  the  equally  spaced  blade  row  of  NB-blades  is  given  by 
Equation  (46) : 
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$ 

s 


1_ 

4tt 


JVS  NB  2 (C  - £.) 

£ 2 (±  

j = l j-1  1 - V 


1 

e 

J 

V 


Qj(n)  dn  - 


- E exp  (;k  (c  - 5 )) 
n=2  J 


F (K  ,v,j) 
sl  on 


on 


Z (K  n) 
o on 


“ 00  F (K  ,v,j) 

- 2 E cos  m(0  - 0.  .)  E exp  (+K  (£  - <;.))  S v" Z (K  pi)  - 

, it  , r mnv  j K m mn 

m=l  J n=l  J mn 

(46) 

Taking  advantage  of  the  characteristics  of  the  equally  spaced  blade  row, 
as  shown  in  Appendix  B,  the  summation  with  respect  to  i is  eliminated 
Also,  the  summation  with  respect  to  m remains  only  when  m is  a multiple 
of  the  blade  number  NB,  then. 


NB  ^ 2CC-5J  1 
*s  = \ C*  - — h / Vn)  dri 

j = l 1 - v v J 


00  F (K  ,v,j) 

1 exp  (;Kon(C  - tj))  ^ Zq 


n=2 


on 


00  00  F (K,  ,v, j) 

- 2 E cos  k (0  - 0 ) E exP(;Kkn(,  - Cj))  ^(K  n)J- 

x=l  J n*l  J kn 

(47) 

Again,  a closed  form  of  the  series  expression  for  the  other  summation 
has  not  been  found  so  far.  The  corresponding  three-dimensional  perturbed 
flow  field  due  to  a series  of  radial  source  lines  of  variable  strength 
replacing  the  equally  spaced  blade  row  of  NB-b lades  is  given  by  the 
following  Equations.  The  upper  and  lower  signs  refer  to  ^ > 0 and 

L - < 0 respectively. 
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u l l k sin  k (9  - 0 ,) 

s9  27TT1  j = l 4=1  13 


JVS  « 


F (K.  ,V,j) 

Z exP  ‘ S3)  “H? Zk  (Kkn^ 


(48) 


n=l 


kn 


“sn=  - JF  ,E  ,<  Z,  K 

' j=l  n=2 


JVS  °° 


F (K  ,v,j)  3Z  (K  n) 
sv  on  o on  ' 


on 


9n 


+ 2 Z cos  k (0-0..)  Z exp  (;K . (C-C ,))  s 

4=1  13  n=l  Kn  3 kn 


Fs  ^Kkn,V'33  9VKfcnn\ 


3n 


(49) 


and 


JVS 


NB 


u » t=-  Z (± 

r /I  tt  ' 


j = l 1-V  V 


E„  exP  ^KonU'Cj3)  Fs(Kon’V>3)  Zo(Konn) 
n=2  J 


± 2 1 cos  k(0-01;j)  Z^exp  ( ( 5 -5^ ) ) Fs(Kkn’V,3)  Zk(Kknn))< 

(50) 

where  u Q,  u , u are  velocity  components  in  the  directions  0,  n,  and 
s0  sf]  s? 

x,  respectively,  and  k = 4 NB, 


3.2  Direct  Problem  in  Axial  Turbomachinery 

The  direct  problem  in  turbomachinery  is  to  determine  the  flow 
field  around  a given  blade  profile.  The  objective  of  this  section  is 
to  provide  a theoretical  formula  for  the  solution  of  the  three- 
dimensional  direct  problem  in  axial  turbomachinery  In  the  distributed 
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singularity  method,  the  strength  of  the  distributed  vortices  and  sources, 
and  the  cascade  mean  velocity  are  not  known,  but  the  configuration  of 
the  blade,  inlet  flow  angle,  stagger  angle  and  spacing  are  given 
Therefore,  the  strength  of  the  singularities  must  be  determined  for  the 
given  blade  geometry,  spacing,  stagger  angle  and  inlet  flow  angle  The 
present  method  is  valid  only  for  a free  vortex  type  blade  row  wi*h 
approximately  radial  elements.  Such  a blade  row  will  be  replaced  by  a 
series  of  radial  vortex  lines  of  constant  strength  and  a series  of 
radial  source  lines  of  variable  strength.  Furthermore,  the  resultant 
flow  field  should  be  irrotational  with  respect  to  the  coordinate  system 
in  which  the  blade  row  is  fixed.  All  the  flow  quantities  are  measured 
with  reference  to  this  coordinate  system.  The  velocities  are  normalized 
by  the  axial  component  of  the  three-dimensional  cascade  mean  velocity 


The  blade  system  can  be  represented  by: 


a)  The  blade  system  has  NB  blades  and  each  blade  is  replaced  by 
JVS  radial  vortex  lines  and  source  lines  respectively. 

b)  Each  vortex  line  is  radial  and  of  constant  strength  along  the 
radial  direction,  and  the  strength  of  the  j-th  vortex  line  for  each 
blade  is  given  by  I\  - 

c)  Each  source  line  is  radial  and  of  variable  strength  along  the 
radial  direction.  The  strength  of  the  j-th  source  line  for  each  blade 
is  represented  by  a p-th  order  polynomial  given. 


% 


(n)  ■ a 


j.o 


A. 

J.1 


(51) 


d)  For  convenience,  the  j-th  vortex  line  and  j-th  source  line  are 


placed  at  the  same  location. 
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Similar  to  the  two-dimensional  solution  of  the  direct  problem  by 

the  distributed  singularity  method,  the  flow  field  due  to  a series  of 

vortex  and  source  lines  must  satisfy  the  kinematic  condition,  the 

closure  condition  for  the  source  distribution,  and  the  Kutta  condition 

First,  a detailed  discussion  of  the  resultant  velocity  and  the  cascade 

# 

mean  velocity  in  the  three-dimensional  flow  is  made  Let  V be  a 

nondimensional  resultant  velocity  vector;  then, 

V = u + u + W , (52) 

vs 

where  u^,  us  are  the  perturbed  velocities  due  to  a series  of  radial 
vortex  lines  and  a series  of  radial  source  lines  respectively,  and  W 
is  the  cascade  mean  velocity.  Let  i,  j,  k be  unit  vectors  in  the 


0,  p,  £ directions  respectively.  Then, 

Jve 


-t  t 

u = l • u Q + i • u + k * u ^ 
v v0  vn  v£ 


(S3) 


and 


u = i,un+j*u  + k • u ^ 
s s0  J sn  s£ 


i • W_  + j • W + k • W„ 

0 J n ? 


(54) 


(55) 


where  uvg,  u^,  u are  the  0,  n»  C components  of  u^.  Similarly,  ucq, 

u , Us^  are  t*ie  r'»  *=>  components  of  ug,  and  Wg,  W^,  W are  the  0, 

—V 

ri,  C components  of  W.  As  described  in  Section  (3.1),  the  cascade  mean 
velocity  in  the  three-dimensional  case  has  only  tangential  and  axial 
components,  and  must  satisfy  the  relation  given  by  Equation  (36),  (37), 
and  (38)  in  Section  (3.1).  Therefore,  the  tangential  component  of  the 
cascade  mean  velocity  at  radius  n is  given  by  the  following  equation: 


W.  = — Wfi 

0 n 0m 


(56) 
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where  = l+v/2,  and  is  the  tang+ntial  component  of  the  cascade 


mean  velocity  at  mid  radius  between  hub  and  outer  casing  wall  The 
axial  component  of  the  cascade  mean  velocity  must  be  constant  every 
where,  and  the  radial  component  must  be  zero  everywhere.  Then, 


W = Wr  =1 
C Cm 


(57) 


and 


W = 0 

n 


(58) 


where  W^m  is  the  axial  component  of  the  cascade  mean  velocity  at  mid 

radius  On  each  cylindrical  surface  of  radius  n,  we  can  define  the 

cascade  mean  flow  angle  as  similar  to  the  cascade  mean  velocity  Let 

fa.)  , (AV„)  and  (a  ) be  respectively  the  cascade  inlet  flow  angle  at 
• l m 0 m mm 

mid  radius,  the  tangential  velocity  change  at  the  exit  at  the  mid 
radius  and  the  cascade  mean  flow  angle  at  the  mid  radius.  Similar  to 
the  two-dimensional  cascade  theory, 


, , "0»  * 
ta"  r 

Cm 


tan  (a  ) + 


(AVm 


mm  W 


(59) 


Cm 


and  also 


2 * 5 * (AVm  = - rt  ' 


(60) 


where  s is  the  spacing  and  T is  the  circulation  along  a circuit  on  the 
cylindrical  surface  at  mid  radius  = l+v/2  (see  Figure  2 and  Figure  4), 


and 


s = 2ttt|  / NB 
m 


JVS 

r = i r. 

t j-i  J 


(61) 


(62) 


r 


k.  i 
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Therefore,  substituting  Equation  (60)  into  (59),  we  get 

tan  (a  ) = tan  (a.)  + I\/2*s*lV  . (63) 

-•  ^ m m i m t Cm  v 1 

In  Equation  (63),  (oc)m  is  given  and  w^m=  1 since  all  the  velocities  are 


normalized  by  this  velocity  component.  Then, K an  (a,  ) is  simply  the 


linear  equation  of  T . From  the  velocity  triangle  (see  Figure  4),  we 
get 


WQ  = W_  • tan  (a  ) 
0m  £m  mm 


tan  (a  ) 
m m 


(64) 


Then,  Equation  (56)  is  rewritten  as 


WQ  = l+v/2n  * tan  (a  ) 

0 v mm 


(65) 


and  also  from  the  definition. 


W = 0 

n 


(66) 


and 


W = Wr  =1  . 


(67) 


Then,  the  components  of  the  cascade  mean  velocity  at  any  radius  are 

represented  in  terms  of  the  cascade  mean  flow  angle  at  mid  radius 

(Equations  (65),  (66),  and  (67)).  Therefore,  the  components  of  the 
. *►  , 

resultant  velocity  vector  V given  in  Equation  (52)  are  rewritten  as 

(68) 


V0  '=  Uv0  + Us0  + 1+V/2n  tan  ^m  * 


V = u + u 
n vn  sn 


(69) 


and 


V = u,+u,+  l 
C vC  sC 


(70) 
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For  representing  the  blade  surface,  it  is  convenient  to  define  two 
planes  such  that  one  is  flat  and  perpendicular  to  the  axial  axis  £ and 
the  other  is  a cylindrical  plane  of  radius  n.  These  are  designated  as 
the  £ -constant -plane  and  the  n -constant -plane  respectively.  As  shown 
in  Figure  5a  and  5b,  let  a be  the  angle  between  the  ^-direction  and  the 
line  AB*  tangential  to  the  blade  surface  at  point  P on  the  n constant 
plane,  and  6 be  the  angle  between  the  radial  line  and  the  tangential 
line  CD  at  a point  P on  the  constant  plane.  The  components  of  the 
unit  normal  vector  n at  a point  P on  the  blade  surface  are  n 0,  n^  and 
n^  in  0,  n and  £ directions  respectively.  Then, 

n = l • nQ  + T • n^  + it  • n^  , (71) 

where 


n = cos  0 cos  a , (72) 


and 


n 


n 


= sin  3 cos  a , 


= cos  3 sin  a 


(73) 


(74) 


Kinematic  Condition: 

The  kinematic  condition  is  that  the  resultant  velocity  vector  must 
be  perpendicular  to  the  vector  n.  Then, 

n • V = 0 . (75) 

Substituting  Equations  (68)  to  (74)  into  (75),  the  kinematic  condition 
is  then  given  by  Equation: 
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cos  S cos  a fu  . + u . + l+v/2n  tan  (a  ) ) 

^ v8  s0  mm 

+■  sin  0 cos  a • (u  + u ) + cos  8 sin  a*(u+u+l)=0, 

v vn  sir  vs  s? 

(76) 


where  a,  Bare  given  by  the  blade  geometry,  and  u^g,  u^,  and  u 
include  to  Tj^g  as  unknown  linear  constants  to  be  determined.  Also, 


u Q,  u , and  u „ incluie  An  , , A.  _,  A,  , 
s0’  sn  s?  1,1  1,2’  1,3 

......  A. 


AAA 
' ’ l,p’  2,1’  2,2’ 

• • • ’ ajvs,p  as 

linear  constants  to  be  determined.  The  term  tan  fa  ) includes  T, , T_, 

^ m m 12 


2,p’ 


’ AJVS,1’  AJVS, 2 ’ 


’ ’ ^JVS  as  l^near  constants  given  in  Equation  (63) 


Closure  Condition: 

Hie  total  source  strength  replacing  the  blade  row  must  be  zero, 


JVS  l JVS  : p 

Z / Q (n)  dn  = Z / ( Z A nS  dp  = 0 (77) 

j=l  V 3 j=l  V i=0 


which  can  be  simplified  to 


JVS  p 

Z Z 
j=l  1=0 


1 


H+l 

v 

1+Jl 


0 


(78) 


Kutta  Condition: 

In  potential  flow  theory,  the  circulation  around  a given  blade 
profile  is  not  unique  without  the  Kutta  condition  even  if  the  kinematic 
condition  is  satisfied.  Therefore,  to  apply  potential  flow  theory  to 
the  actual  flow  around  a blade  profile,  a restriction  on  the  distribution 
of  the  vortices  is  required  based  on  the  Kutta-Joukowski  hypothesis 
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The  Kutta-Joukowski  hypothesis  states  that  the  velocity  must  be  finite 

at  the  trailing  edge.  In  the  distributed  singularity  method,  the 

strength  of  the  vortices  at  the  trailing  edge  must  vanish  to  satisfy 

the  Kutta-Joukowski  condition.  In  the  three-dimensional  direct  problem, 

the  intersections  of  the  blade  and  p-constant  plane  gives  a three- 

dimensional  contour.  However,  the  radial  component  of  the  resultant 

velocity  is  always  normal  to  the  plane  of  this  contour.  Then,  this 

contour  is  specified  by  only  the  axial  and  tangential  components  of  the 

resultant  velocity  and  the  vortex  lines  are  all  radial;  then  the 

circulation  around  these  contours  at  any  radius  is  always  equal  to  T , 

and  the  three-dimensional  Kutta-Joukowski  condition  is  reduced  to  that 

for  each  contour  at  each  radius.  In  Figure  6,  let  points  A and  B be 

■+  •+ 

very  close  to  the  trailing  edge  T on  the  n-constant  plane  V , 

are  the  projections  of  the  resultant  velocity  V on  the  n-constant 

plane  at  points  A and  B respectively.  and  V'A0  are  axial  and 

tangential  components  of  VA.  Similarly,  V and  V are  axial  and 

A DL,  DO 

-► 

tangential  components  of  V . Then,  the  Kutta-Joukowski  condition  is 

D 

->•  -V- 

that  VA  and  should  be  of  the  same  magnitude  and  go  towards  T for  a 
sharp  trailing  edge.  Then, 

VA  cos  e + VAQ  sm  e = cos  uj  + VB0  sin  oj  , (79) 

where  e and  u)  are  given  from  the  blade  geometry.  Equation  (79)  includes 

T and  A . (j  = l to  JVS,  and  SL=  0 to  p)  as  linear  constants.  Now,  the 
1 J > ~ 

kinematic  condition,  closure  condition,  and  the  Kutta-Joukowski 

condition  are  linear  first  order  equations  concerned  with  the  unknown 

constants  T.  and  A . (j  = l to  JVS,  and  £,=0  to  p)  . The  number  of 
] J i 

unknowns  is  JVS  • (p  + 2) - Therefore,  if  JVS  • (p  + 2)  - p blade 


points  are  given  for  the  kinematic  condition,  together  with  one  closure 
condition  and  the  (p  + 1)  Kutta-Joukowski  condition  for  (p  + 1)  n- 
constant -plane,  the  equations  become  simultaneous  first  order  equations 
with  JVS  • (p  + 2) -unknowns.  Hence,  the  strength  of  the  vortex  lines 
and  source  lines  can  be  determined  by  solving  the  simultaneous  equations 
The  velocity  components  are  determined  by  Equations  (52)  to  (55)  using 
Equations  (43)  to  (45)  and  (48)  to  (50).  Also,  the  pressure  coefficient 
are  given  by  the  following  expressions: 


<Vw  = 


lili 

I wl 2 


(80) 


and 


<Vw 

r rr 


mi 

livl2 


(81) 


where  (C  ) and  (C  ) ^ are  the  pressure  coefficients  with  respect  to 
^ ^ m 

the  cascade  mean  velocity  at  radius  q and  the  cascade  mean  velocity  at 


mid  radius  respectively. 
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CHAPTER  IV 

NUMERICAL  TECHNIQUE  AND  SOLUTIONS 


4 , 1 Determining  of  Eigen  Values  and  Related  Functions 

The  three-dimensional  flow  field  due  to  radial  vortex  lines  and 
source  lines  given  in  Equations  (22)  to  (24),  (33)  to  (35),  (43)  to  (45) 
and  (48)  to  (50)  includes  an  infinite  number  of  eigen  values  n and 
related  functions  F (K  ,v) , F (K  ,v,i),  Z (K  q)  and  3 / 3r(  Z (K  q)  , 
These  eigen  values  are  determined  by  Equation  (11),  and  are  the  same 
both  in  the  radial  vortex  system  and  radial  source  line  system  since 
the  boundary  condition  with  respect  to  the  radial  velocity  at  the  hub 
and  the  outer  casing  wall  is  the  same.  However,  the  eigen  values 
related  to  the  0-th  order  Bessel  function  has  meaning  only  for  the 
source  system.  The  first  three  eigen  values  for  each  order  of  Bessel 
function  were  found  step  by  step  and  other  approximate  eigen  values  were 
found  by  extrapolation.  All  eigen  values  are  improved  to  obtain  higher 

_9 

accuracy  by  successive  iteration  up  to  10  , The  eigen  values  were 

obtained  for  Bessel  functions  of  order  m = 0 to  100,  and  the  number  of 
eigen  values  n = 1 to  20  for  each  order  m for  hub-to-tip  ratios  v = 0,3 
and  0,6  (see  Figure  7,  8),  The  eigen  values  in  Reference  (4)  are 
combined  with  the  present  n by  the  following  expression: 

y = v • K 
m,n  m,n 

Now,  the  related  functions  F (K  ,v)  and  F (K  ,v,j)  are  again  given 

v m,n  ' s m,n'  6 6 


by 


1 

/ pZ  (K  n)  dn 
1 m mn 

F (K  , V ) = -4 s 

v mn  fl  2 


I'  qZ  "(K  n)  dq 
1 ^ m mn 


(82) 


the  numerators  of  Equations  (82)  and  (83)  were  resolved  into 
1 

/ dn  (p  = 0 to  3)  and  calculated  by  the  Gaussion  quadra- 

ture formula  to  order  29  for  each  eigen  value  (Reference  5) , which  is  a 
more  accurate  numerical  integral  formula  than  the  trapezoidal  or 
Simpson's  formula  for  the  same  number  of  points.  The  denominators  of 
Equations  (82)  and  (83)  were  calculated  by  integral  formula  for  Bessel 

function  for  each  eigen  value.  Calculation  of  Z (K  q)  and  9/9q 

6 nr  m,n 

Z (K  n)  is  very  time  consuming.  Therefore,  nine  radial  coordinates 
m m,n 

are  selected  from  the  values  calculated  by  the  Gaussian  quadrature 
abscissae  and  q = v,  1.  The  radial  coordinates  for  which  all  related 
functions  and  necessary  Bessel  functions  are  prepared  are  given  in 
Tables  1,  2,  and  3- 

Table  1,  Radial  Coordinates  Used  for  Computation. 


k 1 2 

3 

4 

5 

6 

7 8 9 

nk  v f(x?) 

fUi0) 

fCx^) 

(l+v)/2 

£<Xi7) 

f(x2Q)  f(x23)  1.0 

where 

f(x^)  = ((l+v+(l-v)  • xi)/2.0  (84) 

and 
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x?  = - x23  = -0.7524629 

x10  = - x20  = -0.5075930  (85) 

X13  = ” xi7  = -0.2113523 

These  correspond  to  the  i-th  Gaussion  abscissa  of  order  29.  The  values 
for  v = 0.6  and  0.3  are  tabulated  in  Table  2 and  3 respectively. 


Table  2.  Radial  Coordinates  Used  for  Computation  (v  = 0.6), 


k 

1 

2 

3 

4 

5 

6 

7 

8 

9 

nk 

0,60000 

0,64951 

0.69848 

0,75773 

0.80000 

0.84227 

0.90152 

0.95049 

1.00000 

Table  3.  Radial  Coordinates  Used  for  Computation  (v  = 0.3). 


k 

1 

2 

3 

4 

5 

6 

7 

8 

9 

nk 

0.30000 

0.38664 

0.47234 

0.57603 

0.65000 

0,72397 

0,82766 

0.91336 

1.00000 

The  values  Z (K  r|)  for  the  above  mentioned  radial  coordinates  are 
m m,n 

already  calculated  in  the  process  of  evaluating  the  numerical  integral 
1 

/ nPZ  (K  n)  dn  and  also  the  values  of  9/3n  Z (K  n)  were  calculated 
' m m,n  m m,n 

by  known  values  using  the  recurrence  formula  for  the  Bessel  function. 

To  avoid  duplicate  calculation  of  the  Bessel  function  in  the  main  program, 
all  eigen  values  and  related  functions  were  stored  on  a magnetic  tape 
for  m = 0 to  100,  n = 1 to  20,  v = 0,3  and  0,6  for  the  radial 
coordinates  given  in  Tables  2 and  3,  Then,  the  computer  program  for 
the  three-dimensional  potential  flow  field  due  to  an  arbitrary 
series  of  radial  vortejc  lines  of  constant  strength  and  a series  of 
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radial  source  lines  of  variable  strength  can  be  used  for  hub-to-tip 
ratios  v = 0,3  and  0.6  (see  Reference  6).  Furthermore,  if  the  eigen 
values  and  related  functions  are  obtained  in  a similar  manner  for  any 
other  hub-to-tip  ratio,  the  above  computer  program  is  valid  for  that 
hub-to-tip  ratio  as  well.  The  computer  programs  for  the  determination 
of  the  eigen  values  and  related  Bessel  functions  are  given  in  Reference 
6. 


4 .2  Numerical  Calculation  of  the  Three- 

Dimensional  Perturbed  Flow  Field 

4.2.1  Convergency  of  the  Solution  of 
the  Three-Dimensional  Perturbed 
Flow  Field 

The  three-dimensional  perturbed  flow  field  due  to  a series  of 
radial  vortex  lines  of  constant  strength  and  a series  of  radial  source 
lines  of  variable  strength  given  by  Equations  (43)  to  (45)  and  (48)  to 
(50)  include  an  infinite  number  of  summations  with  respect  to  m and  n. 
In  order  to  examine  the  convergency  of  the  series  expressions,  a single 
radial  vortex  line  of  strength  unity  located  at  £ = 0,  0 = 0°  and  a 
single  radial  source  line  of  strength  unity  located  at  the  same  point 
were  examined  for  v = 0.6  separately.  The  contributions  of  m and  n to 
the  convergency  of  the  tangential  velocity  at  £ = 0,05,  0.1,  0-2  and 
0 = 30°  for  n * 08  due  to  a single  radial  vortex  line  are  shown  in 
Figure  9 The  corresponding  graphs  for  a single  radial  source  line 
are  plotted  in  Figure  10,  It  is  clear  that  the  convergency  of  the 
series  expression  strongly  depends  on  m,  not  on  n.  For  hub-to-tip 
ratio  v = 06,  the  required  order  of  Bessel  function  m is  approximately 
m > 60  for  ls|  < 01,  m > 20  for  |^|  >0.2  and  n = 6 is  enough  in  any 
case-  Therefore,  these  results  suggest  that  the  summation  based  on 
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only  n = 1 for  each  m will  give  a good  result  for  large  |?|  if  m is 
taken  large  enough  However,  the  convergency  becomes  poor  at  points 
close  to  the  singular  point  | ^ | <0 - 05  for  v = 0,6  even  if  m is  taken  up 
to  100  This  situation  is  more  severe  in  the  case  of  smaller  hub-to- 
tip  ratio.  It  seems  the  smaller  the  hub-to-tip  ratio  is,  the  more  the 
eigen  values  are  required.  Mathematically,  the  three-dimensional 
tangential  velocity  components  due  to  a vortex  line  and  the  axial 
velocity  component  due  to  a source  line  located  at  £ = 0 and  0=0° 
must  vanish  at  £ = 0 for  any  0 except  0=0°  for  any  hub-to-tip  ratio 
Therefore,  taking  advantage  of  these  characteristics  of  the  series 
expression,  some  improvement  to  avoid  the  above  difficulties  can  be 
made  by  further  analytical  considerations  described  in  Appendix  D, 

This  approach  will  be  pursued  later. 

422  Three  Dimensional  Perturbed  Flow 

Field  Due  to  a Single  Radial  Vortex 
Line  of  Strength  Unity  and  Its 
Comparison  with  the  Two-Dimensional 
Solution 

The  three-dimensional  perturbed  velocity  components  due  to  a single 
radial  vortex  line  of  strength  unity  located  at  = 0 and  0^  ^ = 0° 
were  calculated  by  Equations  (43)  to  (45)  for  hub-to-tip  ratios  V = 0 6 
and  03,  The  tangential  velocity  component  at  the  hub  due  to  a single 
radial  vortex  line  is  compared  with  the  two-dimensional  values  in 
Figure  11  for  v = 0 6 and  Figure  25  for  v = 0.3,  The  values  in  the  poor 
convergency  region  have  been  excluded.  The  discrepancy  between  the 
three-dimensional  and  two-dimensional  solutions  is  large  for  small  0 
and  the  maximum  difference  is  about  25%  of  two-dimensional  values  for 
v = 0 6 The  qualitative  characteristics  of  the  three-dimensional 
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tangential  velocity  component  are  very  similar  to  those  of  two- 
dimensional  values,  and  at  £ = 0,  X,  -*■  ± °°,  the  velocities  have  the 
same  limiting  values  mathematically  (see  Appendix  D) . The  effects  of 
the  hub-to-tip  ratio  v on  the  tangential  velocity  component  is 
significant  and  the  maximum  difference  from  the  two-dimensional  solution 
is  about  30%  in  the  small  0 region  for  v = 0,3  (see  Figure  25).  However, 
accurate  three-dimensional  values  at  points  close  to  the  singular  points 
could  not  be  obtained  due  to  the  poor  convergency  of  the  series 
expression  as  mentioned  in  Section  4,2.1. 

The  tangential  velocity  components  at  the  mid  radius  and  tip  due 
to  a single  radial  vortex  line  of  strength  unity  are  shown  in  Figures 
12  and  13  for  hub-to-tip  ratio  v = 0.6,  and  corresponding  graphs  for 
v = 0 3 are  shown  in  Figures  26  and  27.  At  mid  radius  (Figure  12), 
the  discrepancies  between  the  three-dimensional  and  two-dimensional 
solutions  vanish  when  the  hub-to-tip  ratio  v = 06,  but  very  small 
differences  were  found  for  v = 03  (see  Figure  26).  For  both  v = 0 6 
and  03,  the  three-dimensional  tangential  velocity  components  at  the 
hub  foi  small  0 are  greater  than  the  two-dimensional  values,  while  those 
at  the  tip  for  small  0 are  smaller  than  the  two-dimensional  values  (in 
absolute  values)  These  three-dimensional  inviscid  effects  may 
partially  account  for  the  higher  (than  the  two-dimensional  values)  lift 
coefficients  measured  at  the  root  and  lower  lift  coefficients  measured 
at  the  tip  of  a pump  rotor  by  Oshima  et  al,  (Reference  7) 

The  axial  velocity  components  at  the  hub,  mid  radius  and  tip  due  to 
a single  radial  vortex  line  of  unit  strength  are  shown  in  Figures  14, 

15,  and  16  respectively,  for  hub-to-tip  ratio  v = 0.6  Corresponding 
graphs  for  v = 0,3  are  shown  in  Figures  28,  29,  and  30  respectively. 
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Similar  to  the  tangential  velocity  components,  the  discrepancy  between 
the  three-dimensional  and  two-dimensional  axial  velocity  components 
vanish  at  mid  radius  for  V = 0,6  (see  Figure  15),  and  very  small 
differences  were  found  at  mid  radius  for  v = 0-3  (see  Figure  29)-  Also, 
the  three-dimensional  axial  velocity  components  at  the  hub  for  small  £ 

(<  0 2)  are  smaller  than  two-dimensional  values  (see  Figures  14  and  28) 
for  both  hub-to-tip  ratio  v = 0.6  and  03.  The  magnitude  of  the 
discrepancy  between  the  three-dimensional  and  two-dimensional  axial 
velocity  components  for  small  £ region  are  10  to  15%  for  v = 0,6. 

However,  for  v = 0,3,  for  the  small  £ region,  the  numerical  results  of 
the  three-dimensional  axial  velocity  components  are  poor.  Thus,  the 
discrepancy  is  not  clear,  but  a larger  discrepancy  than  that  for  V = 06 
may  be  expected 

The  radial  velocity  components  at  the  mid  radius  due  to  a single 
radial  vortex  line  of  strength  unity  are  shown  in  Figure  17  for  v = 0.6 
in  Figure  31  for  v = 0-3-  The  radial  velocity  components  vanish  at  the 
hub  and  the  outer  casing  wall  The  typical  radial  velocity  distributions 
are  shown  at  the  mid  radius-  However,  the  mid  radius  is  not  necessarily 
the  point  of  maximum  radial  velocity.  The  maximum  radial  velocity 
usually  occurs  at  a smaller  radial  coordinate  than  mid  radius.  From  the 
comparison  of  the  radial  velocity  components  at  mid  radius  for  v = 0 6 
and  v = 03,  the  value  for  v = 0,3  is  about  3 or  4 times  larger  than 
that  for  v = 0,6.  From  analytical  considerations,  it  can  be  shown  that 
the  radial  velocity  components  must  vanish  at  X,  = 0 for  all  0 except 
0=0°.  The  graphs  do  not  show  the  small  C region  due  to  the  lack  of 


convergency  of  the  series  expression. 
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It  is  clear  from  the  plots  shown  in  Figures  11  to  17  and  Figures 
25  to  31  that  the  effects  of  the  three-dimensionality  are  amplified  at 
small  hub-to-tip  ratio  (v) , The  order  of  the  magnitude  of  the  radial 
velocity  components  is  almost  10%  of  the  tangential  and  axial  velocity 
components  for  hub-to-tip  ratio  v = 0.6,  but  for  v = 0.3,  it  is  more 
than  15%.  With  regard  to  tangential  and  axial  velocity  components,  the 
discrepancy  between  two-dimensional  and  three-dimensional  solutions  is 
opposite  in  relative  magnitude  at  the  hub  and  tip  locations  for  any  hub- 
t.o-tip  ratio  This  discrepancy  increases  with  decrease  in  hub-to-tip 
ratio  Finally,  the  three-dimensional  numerical  results  were  compared 
with  References  (2)  and  (4),  The  present  numerical  results  for  v = 06 
are  coincident  with  the  results  of  Reference  (4) , but  different  from  the 
results  of  Reference  (2) , However,  no  other  numerical  results  are 
available  for  a small  hub-to-tip  ratio  such  as  v = 0-3  for  comparison 

423  Three-Dimensional  Perturbed  Flow  Field 
Due  to  a Single  Radial  Source  Line  of 
Strength  Unity  and  Its  Comparison  with 
the  Two-Dimensional  Solution 

Numerical  calculation  of  the  three-dimensional  flow  field  due  to  a 
radial  source  line  has  not  been  reported  earlier-  Even  though  the  flow 
field  around  the  blades  depends  mainly  on  the  vortex  distribution,  the 
three-dimensional  effects  due  to  a source  line  (representing  the  thick- 
ness) can  be  of  the  same  order  as  the  three-dimensional  effects  due  to 
a vortex  line  This  is  especially  true  when  there  is  appreciable 
radial  variation  in  blade  thickness.  The  three-dimensional  perturbed 
velocity  components  due  to  a single  radial  source  line  of  strength 
unity  located  at  = 0 and  0^  ^ = 0°  were  calculated  from  Equations 

(48),  (49)  and  (50)  for  the  hub-to-tip  ratios  v = 0.6  and  0,3.  The 


computer  program  can  be  used  for  the  determination  of  the  three- 
dimensional  flow  field  due  to  a series  of  radial  source  lines  of 
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variable  strength  represented  by  the  third -order  polynomials,  but  only 
the  case  of  constant  strength  (=  1.0)  is  shown  in  Figures  18  to  24  and 
Figures  32  to  38  together  with  two-dimensional  solutions.  The  effects 
of  the  variable  strength  of  the  radial  source  lines  will  be  shown  in 
Section  (4.3.4) . 

The  tangential  velocity  components  at  the  hub,  mid  radius  and  tip 
due  to  a single  radial  source  line  of  strength  unity  are  shown  in 
Figures  18,  19,  and  20  for  a hub-to-tip  ratio  v = 06,  and  Figures  32, 
33,  and  34  for  V = 0.3.  Consider  the  two-dimensional  flow  field  due  to 
a vortex  of  strength  unity  and  a source  of  strength  unity.  The  magni- 
tude of  the  two-dimensional  tangential  velocity  components  due  to  a 
vortex  is  the  same  as  the  two-dimensional  axial  velocity  components  due 
to  a source  except  for  the  difference  of  sign-  Similarly,  the  two- 
dimensional  axial  velocity  components  due  to  a vortex  are  the  same 
magnitude  as  the  two-dimensional  tangential  velocity  components  due  to 
a source  The  relative  characteristics  of  the  three-dimensional 
perturbed  flow  fields  due  to  a single  radial  vortex  line  of  strength 
unity  and  those  due  to  a single  radial  source  line  of  strength  unity  are 
very  similar  to  the  relative  characteristics  of  the  two-dimensional  case 
as  mentioned  above.  Hence,  the  three-dimensional  tangential  velocity 


components  at  the  hub  due  to  a single  radial  source  line  (see  Figure 
18)  are  quite  similar  to  the  axial  velocity  components  at  the  hub  due 
to  a single  radial  vortex  line  (see  Figure  14) . However,  the  discrepancy 
between  two-dimensional  and  three-dimensional  solutions  are  opposite  in 


38 


A 


line  is  larger  than  the  two-dimensional  value  (see  Figure  18),  but  the 
axial  velocity  at  the  hub  due  to  a vortex  line  is  smaller  than  the  two- 
dimensional  value  (see  Figure  14) . Similarly,  the  tangential  velocity 
components  at  the  mid  radius  (see  Figure  19)  and  the  tip  (see  Figure  20) 
due  to  a source  line  are  very  similar  to  the  axial  velocity  components 
at  mid  radius  (see  Figure  15)  and  the  tip  (see  Figure  16)  due  to  a 
vortex  line.  It  is  clear  from  Figure  19  that  the  tangential  velocity 
components  at  the  mid  radius  due  to  a source  line  are  the  same  as  the 
two-dimensional  values  as  in  the  case  of  a vortex  (see  Figure  15) . The 
tangential  velocity  components  at  the  tip  due  to  a source  line  shown  in 
Figure  20  are  smaller  than  the  two-dimensional  values;  however,  in  the 
corresponding  graph  of  Figure  16,  the  axial  velocity  due  to  a vortex 
line  is  larger  than  that  of  the  two-dimensional  value  Comparing  the 
tangential  velocity  components  at  the  hub  and  at  the  tip  due  to  a source 
line  we  find  that  at  the  hub  the  three-dimensional  values  are  larger 
than  the  corresponding  two-dimensional  Values  while  at  the  tip  they  are 
smaller  than  the  two-dimensional  values  The  above  discussion  for  the 
three-dimensional  tangential  velocity  components  due  to  a source  line  for 
the  hub-to-tip  ratio  v = 06  is  qualitatively  valid  for  the  case  of  the 
hub-to-tip  ratio  v = 0-3.  Specifically,  Figures  32,  33,  and  34  are 
similar  to  Figures  28,  29,  and  30  respectively,  and  the  characteristics 
of  the  three-dimensional  tangential  velocity  components  relative  to  the 
two-dimensional  values  for  v = 0.3  due  to  a source  line  are  the  same 
as  those  for  hub-to-tip  ratio  v = 0,6.  The  relative  behavior  of  the 
tangential  velocity  for  v = 0.3  due  to  a source  line  at  the  hub  and  tip 
is  also  quite  similar  to  that  for  v = 0,6,  However,  the  magnitude  of 
the  discrepancy  between  the  three-dimensional  tangential  velocity 
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components  due  to  a source  line  and  those  for  the  two-dimensional  case 
is  amplified  at  small  hub-to-tip  ratio  (v) . At  mid  radius,  the 
tangential  velocity  components  for  v = 0,3  due  to  a source  line  are 
slightly  different  from  the  two-dimensional  values 

The  axial  velocity  components  at  the  hub,  mid  radius  and  tip  due  to 
a single  radial  source  line  of  strength  unity  are  shown  in  Figures  21, 

22,  and  23  for  hub-to-tip  ratio  V = 0,6,  and  in  Figures  35,  36,  and  37 
for  v = 0,3  respectively.  The  axial  velocity  at  the  hub  due  to  a source 
line  for  v = 0,6  shown  in  Figure  21  is  similar  to  the  tangential  velocity 
at  the  hub  due  to  a vortex  line  except  for  the  difference  of  sign.  In 
the  same  way.  Figures  22  and  23  are  similar  to  Figures  12  and  13 
respectively.  The  axial  velocity  at  mid  radius  due  to  a source  line 
shown  in  Figure  22  is  the  same  as  the  tangential  velocity  due  to  a vortex 
line,  and  the  discrepancies  between  the  three-dimensional  and  two- 
dimensional  values  vanish.  The  relative  behavior  of  the  axial  velocity 
due  to  a source  line  at  the  hub  and  tip  is  such  that,  in  the  small  £ 
region,  the  three-dimensional  axial  velocity  components  at  the  hub  are 
smaller  than  those  for  the  two-dimensional  case,  but  at  the  tip  the  three- 
dimensional  values  are  larger  than  the  two-dimensional  values  (see  Figures 
21  and  23)  For  the  case  of  the  small  hub-to-tip  ratio  v = 0.3,  the 
earlier  comments  are  still  valid,  but  the  discrepancies  between  the  three- 
dimensional  axial  velocity  components  due  to  a source  line  and  the  corres- 
ponding two-dimensional  values  are  amplified  by  the  small  hub-to-tip 
ratio  (see  Figures  35,  36,  and  37).  Then,  axial  velocity  at  the  mid 
radius  for  v = 03  due  to  a source  line  is  slightly  different  from  the 
two -dimensional  axial  velocity  component.  This  may  be  partly  due  to  poor 
convergency  of  series  expressions,  but  the  remaining  part  of  the 
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discrepancies  may  be  due  to  essential  characteristics  of  small  hub-to- 
tip  ratio  since  mid  radius  is  not  necessarily  the  point  of  the  same 
velocity  as  two-dimensional  theory.  The  closed  form  of  the  tangential 
and  radial  velocity  component  due  to  a single  radial  source  line  of 
strength  unity  is  not  known  so  far,  but  the  axial  velocity  component  due 
to  a source  line  vanishes  at  £ = 0 for  any  9 except  9 = 0°,  similar  to 
the  two-dimensional  analysis.  However,  the  limiting  value  at  is 

different  from  the  two-dimensional  solution  since  at  (u  ) ^=±1/ 

C2rt(l+V) ) in  the  three-dimensional  case  and  (us^)2^=±l/(4Trn)  in  two- 
dimensional  flow.  Therefore,  these  limiting  values  are  identical  only 
for  the  mid  radius  (n  = (l+v)/2),  as  shown  in  Figures  35,  36,  and  37, 

The  radial  velocity  components  at  the  typical  radial  coordinate 
r - (l*v)/2  due  to  a single  ladial  source  line  of  strength  unity  are 
shown  in  Figure  24  for  v = 0.6  and  in  Figure  38  for  V = 0,3,  respectively 
The  order  of  magnitude  of  the  radial  velocity  components  is  about  10%  to 
15%  cf  the  corresponding  three-dimensional  tangential  and  axial  velocity 
components  for  v = 0,6  and  0,3,  respectively.  The  behavior  of  the  radial 
velocity  components  due  to  a source  line  are  completely  different  from 
those  due  to  a vortex  line, 

4 3 Numerical  Solution,  of  Turbomachinery  Flow 

4 3,1  Lifting  Line  Model 

The  purpose  of  ^his  section  is  to  merely  point  out  the  extent  of 
three-dimensional  effects  due  to  hub-to-tip  ratio  and  the  number  of  the 
blades  NB 


NB  radial  vortex  lines  of  strength  unity  are  located  at  £ = 0-plane 
spaced  equally  along  the  circumference.  The  coordinate  of  the  i-th 


41 


radial  vortex  line  is  given  by  X,  - 0,  0 = 360°/NB  (i-1)  . The  calculated 

tangential  velocity  distributions  are  given  in  Figures  39  and  40,  for 

8=2°  and  9°  representing  regions  very  close  to  and  far  away  from  the 

vortex  lines  These  two  angular  locations  are  smaller  than  the  angle  of 

the  one  spacing  for  any  NB,  but  in  the  case  of  NB  = 20,  9=9°  represents 

the.  mid  passage  The  tangential  velocity  u^  has  the  maximum  absolute 

3d 

value  with  the  maximum  difference  from  two-dimensional  solution  at  0 = 0. 

In  Figure  39,  the  ratio  (uvg  /uvq  ) •’LS  plotted  along  axial  axis-£.  It 

2d  3d 

is  evident  that  the  three-dimensional  effects  are  larger  for  a smaller 

number  of  blades  and  the  ratios  of  the  two  solutions  are  almost  unity 

for  NB  = 20,  In  other  words,  the  value  of  (u  Q ) will  be  closer  to  the 

V03d 

twc -dimensional  solution  for  a larger  number  of  the  blades.  The  ratio 

(u  /u  ^ ) will  approach  unity  far  downstream  and  far  upstream  from 

2d  3d 

the  vortex  lines.  The  effects  of  the  small  hub-to-tip  ratio  (v  = 0,3) 

on  (u^g  /Up  ) is  large  compared  with  those  for  the  large  hub-to-tip 
2d  3d 

ratio  (v  = 0,6),  The  two-dimensional  solution  gives  an  underestimated 

absolute  value  of  the  tangential  velocity  at  the  hub  and  an  overestimated 

value  at  the  tip  than  the  three-dimensional  solution,  and  this  is  true 

for  any  number  of  the  blades  NB  and  for  all  hub-to-tip  ratios.  However, 

at  mid  radius  (n  = l+v/2) , the  ratio  (u^g  /u  g ) is  very  close  to  unity 

2d  3d 

for  all  the  cases.  The  maximum  discrepancy  between  the  two  solutions 
is  as  much  as  20%  to  30%,  especially  for  low  hub-to-tip  ratio  and  large 
blade  spacing.  It  is  evident  that  two-dimensional  theories  are  inac- 
curate for  low  solidity  and  low  hub-to-tip  ratio  blading  used  in  certain 

applications.  In  Figure  41,  the  behavior  of  the  ratio  (uyg  /uvg  ) for 

2d  3d 

NB  - 3 at  1/2,  1/4,  1/8,  and  1/16  of  the  blade  spacing,  which  corresponds 

to  0 = 60°,  30°,  15°,  and  7.5°  respectively,  is  shown.  The  largest 

departure  of  the  ratio  (u  . /u  a ) from  unity  is  at  mid  passage 

V02d  V03d 
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(8  = 60°},  On  the  other  hand,  the  largest  difference  of  |u  „ - u _ I 

V02d  V03d 

is  1/16  spacing  point  because  of  the  short  distance  from  the  vortex  lines 
It  should  be  pointed  out  that  the  distribution  of  the  axial 
velocity  component  due  to  equally  spaced  vortex  lines  on  the  £ = 0-plane 
decreases  rapidly  to  zero  for  both  two-  and  three-dimensional  solutions, 
with  respect  to  the  increase  in  £,  This  is  especially  true  for  NB  = 20, 
Also,  the  order  of  magnitude  of  the  induced  axial  velocity  due  to 
vortex  lines  is  usually  smaller  than  the  corresponding  tangential 
velocity,  especially  for  NB  = 20,  The  axial  velocity  is  almost  zero  at 
the  hub  for  small  hub-to-tip  ratios  everywhere.  The  simple  comparison 
of  the  ratio  fu  „ /u  „ ) when  the  absolute  value  of  the  axial  velocity 

v^2d  VSd 

is  small  leads  to  misleading  conclusions  on  the  effects  of  NB.  The 
comparison  of  the  effects  of  NB  were  made  at  the  aoint  £ = 0.1,  where 
the  absolute  value  of  axial  velocity  is  relatively  large.  The  calcula- 
ted axial  velocity  ratios  (u  , /u  ) on  the  plane  of  £ = 0,1  are 

v^2d  vhd 

shown  in  Figures  42a,  42b,  and  42c  for  V = 0,6  and  V = 0,3  at  the  hub, 
mid  radius  and  tip,  respectively,  where  s is  the  spacing  angle  and  is 
equal  to  360°/NB.  The  discrepancy  between  the  two  solutions  is  larger 
for  small  number  of  the  blades  at  the  hub  for  both  v = 0,6  and  0,3, 

The  behavior  of  the  ratio  (u  „ /u  „ ) at  the  tip  seems  to  be  different 

v^2d  vSd 

from  that  at  the  hub;  however,  the  order  of  magnitude  of  the  axial 
velocity  at  the  tip  for  NB  = 20  is  about  1/10  of  the  corresponding 
tangential  velocity  component.  On  the  other  hand,  the  order  of  the 
magnitude  of  the  axial  velocity  at  the  tip  for  NB  = 1 is  almost  the 
same  as  that  of  the  corresponding  tangential  velocity.  Therefore, 
three-dimensional  effects  are  large  for  small  hub-to-tip  ratios  with 
fewer  blades  The  two-dimensional  solution  gives  overestimated  absolute 


w 


values  at  the  hub  and  underestimated  values  at  the  tip  for  small  number 

of  the  blades.  The  ratio  (u  /u  ) at  the  mid  radius  is  close  to 

VQ2d  V^3d 

unity  and  the  effects  of  NB  and  hub-to-tip  ratio  v are  relatively 
smaller  at  the  hub  and  tip. 

The  radial  velocity  distributions  are  shown  in  Figure  43  for  9=2° 
and  in  Figure  44  for  0 = 9°  at  mid  radius.  The  radial  velocity  is 
normaliz  J with  respect  to  the  two-dimensional  axial  velocity  at  corres- 
ponding points  since  there  is  no  radial  velocity  component  in  a two- 
dimensional  flow.  Comparing  Figures  43  and  44,  we  find  the  magnitude 
of  the  radial  velocity  ratio  at  two  different  angular  locations  very 
close.  The  hub-to-tip  ratio  has  a dominant  effect  compared  with  the 
number  of  the  blades.  Furthermore  the  radial  velocities  are  about  10% 
to  40%  of  the  two-dimensional  axial  velocity  indicating  the  extent  of 
three-dimensionality  in  the  flow. 


4, 3- 2 Source  Line  Model 

The  induced  velocity  due  to  NB  radial  source  lines  of  unit 

strength  spaced  equally  at  ? = 0-plane  is  to  be  reported.  The  objective 

of  this  investigation  is  to  examine  the  effects  of  the  number  of  the 

blades  NB  and  the  hub-to-tip  ratio  v on  the  flow  field  induced  by 

equally  spaced  NB-radial  source  lines  simulating  the  blade  thickness. 

Generally,  the  flow  field  induced  by  source  lines  are  orthogonal 

relative  to  those  induced  by  vortex  lines  of  the  same  magnitude. 

Namely,  the  distribution  of  the  tangential  velocity  u and  axial 

S03d 

velocity  u due  to  source  lines  are  similar  to  the  distribution  of 


the  axial  velocity  u and  tangential  velocity  u Q due  to  vortex 

VSd  V03d 

lines  respectively 
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The  effects  of  the  NB  on  the  tangential  velocity  ratio  (u  . /u  . ) 

S02d  S03d 

aie  plotted  for  NB  = 1,  3,  5 and  20  at  £ = 0.1  in  Figures  45a,  45b,  and 

45c,  The  order  of  the  magnitude  of  the  tangential  velocity  due  to 

source  lines  u „ is  usually  smaller  than  the  corresponding  axial 
S03d 

velocity  due  to  source  lines.  This  is  especially  true  for  NB  = 20, 

where  tangential  velocity  ugg  is  almost  zero.  Therefore,  from 

Figures  45a,  45b,  and  45c,  the  effects  of  the  NB  on  the  ratio  (u  , / 

S02d 

u Q ) must  be  carefully  considered  since  for  extremely  large  NB,  the 
S03d 

absolute  value  of  the  tangential  velocity  is  very  small  while  the  ratio 
difference  is  quite  large.  Alternatively,  for  smaller  NB,  the  absolute 
value  is  relatively  larger  and  the  ratio  difference  is  relatively  smaller 
It  is  clear  that  the  effects  of  the  hub-to-tip  ratio  are  dominant, 
unlike  the  effect  of  NB,  on  the  three-dimensionality  of  the  flow  field 
As  mentioned  in  the  discussion  of  axial  velocity  induced  by  a 
single  radial  source  line  (Figure  35)  spanning  an  annulus,  the  three- 
dimensional  solution  is  very  different  from  the  two-dimensional  solution 
at  far  up  and  down  stream.  The  two-dimensional  solution  is  the  solution 
expressed  m a cylindrical  coordinate  system-  The  two-dimensional 
unrolled  solution  has  no  interference  effects  due  to  the  flow  field  at 


other  radii  (r,)  , At  far  up  and  down  stream  (£  -*•  ± <*>),  the  ratio 
u /u  , is  independent  of  NB.  In  the  case  of  the  tangential 

S?2d  5Sd 

\elocity  u n and  axial  velocity  u , induced  by  vortex  lines,  and 
V03d  V?3d 

the  tangential  velocity  u „ induced  by  source  lines,  the  ratio  of 

S03d 

these  values  to  the  two-dimensional  solution  is  unity  at  far  up  and 


down  stream  (£  -+  ± °°)  since  u 


and  u 


vanish  at  C = - 00  as  the 


two-dimensional  solution  does.  By  the  irrotationality  condition,  fluid 


w 


layer  at  any  radius  r)  has  the  same  amount  of  circulation  given  by 
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NB*r/2  at  £ = ± 00 . However,  in  the  case  of  the  axial  velocity  u , 

SSd 

induced  by  source  lines,  the  fluid  layer  at  any  radius  n has  an  identi- 
cal magnitude  given  by 


(«  r ) 


2tt(1+v) 


But  the  two-dimensional  value  u „ is  a function  of  the  radius  n even 

S4d 

at  £ = ±°°  since. 


(Us?  3 

s4d 


NB  • Q 
4 v r) 


Hence,  the  ratio  (u  , /u  „ ) has  the  value 

S^2d  SSd 


(US5  /Us£  3 
S^2d  S^3d 


1 + v 
2n 


The  effects  of  the  NB  on  the  axial  velocity  ratio  (u  „ /u  ) are 

sC2d  S?3d 

shown  in  Figure  46  along  the  £ axis  for  several  NB  at  0 = 2 , a region 

relatively  near  the  singular  points.  The  effect  of  the  number  of  the 

blades  is  large  for  large  NB  near  the  radial  source  lines,  but  at  the 

far  down  stream,  the  ratio  is  independent  of  NB  and  has  the  value  1+v/ 

2r\  as  mentioned  above.  The  order  of  the  magnitude  of  the  axial  velocity 

u induced  by  source  lines  is  usually  larger  than  that  of  the  corres- 

S '3d 

ponding  tangential  velocity  u „ which  goes  to  zero  rapidly,  and  this 

Sy3d 

is  also  true  for  the  two-dimensional  case.  The  difference  between  the 
two  solutions  is  great  for  small  hub-to-tip  ratio  (v  = 0,3),  especially 
at  the  hub  However,  at  mid  radius,  the  difference  between  two 


solutions  is  very  small  compared  with  the  differences  at  the  hub  and 
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Finally,  the  radial  velocity  distribution  at  the  mid  radius  induced 

by  NB  radial  source  line  of  unit  strength  located  at  the  plane  £ = 0 is 

calculated  for  V = 0.6  and  0.3  at  9 = 2°  and  is  shown  in  Figure  47  The 

radial  velocity  distributions  are  plotted  in  terms  of  the  ratio  radial 

velocity  to  two-dimensional  axial  velocity  component  (u  /u  ) . As 

sn  s c.2d 

indicated  earlier,  the  radial  velocity  due  to  vortex  lines  does  not 
change  its  sign  with  respect  to  ? or  9 . On  the  other  hand,  the  radial 
velocity  due  to  source  lines  changes  its  sign  with  respect  to  C or  9 
The  plot  shown  in  Figure  47  shows  the  effects  of  the  hub-to-tip  ratio  and 
the  number  of  the  blades  on  the  radial  velocity  induced  in  three- 
dimensional  flow. 


4.3,3  Lifting  Surface  Model  Including 
Thickness  Effects 

In  order  to  examine  the  effects  of  the  stagger  angle,  number  of  the 
blades,  thickness  and  hub-to-tip  ratio,  a system  of  NB  series  of  radial 
vertex  and  source  lines  are  considered  for  several  simple  models  The 
model  consists  of  NB-equally  spaced  blades  and  each  blade  is  represented 
by  eight  vortex  and  source  lines  arranged  equidistant  along  the  camber 
lire  of  the  blades  (see  Figures  48  and  49).  In  the  present  calculation, 
the  camber  is  assumed  to  be  zero  and  the  calculations  are  carried  out 
fer  s*-agg°r  angles  (at  mid  radius)  of  0°,  18°,  36°,  and  54°.  The  nondi- 


mentional  chord  length  of  the  blades  C is  equal  to  2.0  (twice  the  outer 
casing  radius)  on  the  cylindrical  plane  of  radius  l+v/2  (see  Figure  49) 


The  distribution  of  the  strength  of  series  vortex  and  source  lines  Ti , 

Q^,  j = 1 , 2 , ,8  representing  the  blades  are  assumed  to  be  the 

came  at  all  stagger  angles  A,  number  of  the  blades  NB  and  hub-to-tip 


ratio  v Even  though  this  is  not  true  for  a practical  blade 
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configuration,  the  main  objective  of  this  investigation  is  to  understand 

the  extent  of  three-dimensional  effects  and,  from  this  point  of  view, 

these  calculations  are  very  useful.  The  values  of  F,  and  Q.  used  in  this 

JVS  3 3 

calculation  are  given  in  Table  4,  where  Z is  equal  to  -1.0  according  to 
JVS  j=l 

our  definition  and^Z^  is  zero,  satisfying  the  closure  condition. 

Table  4.  Strength  of  Vortex  and  Source  Configuration  Used. 


j-1 

3=2 

j=3 

j=4 

j=5 

j = 6 

j = 7 

3=8 

r.  -o  44000 
j 

-0,16843 

-0,11862 

-0,09069 

-0,07053 

-0. 05360 

-0,03824 

-0,01987 

Q 0.44559 

1 

0. 05441 

-0  03425 

-0,08085 

-0.10545 

-0,11279 

-0,10301 

-0.06366 

The  location  of  the  vortex  and  source  lines  are  specified  only  by  the 
tangential  and  axial  coordinate  (0,  £)  since  the  vortex  and  source  lines 
are  assumed  to  be  radial.  The  stagger  angle  A,  the  coordinate  of  the 
vortex  and  source  lines  and  the  calculation  points  in  (£'  -plane 

and  those  in  (£', 9) -plane  are  shown  in  Figure  49  and  in  Figure  50,  where 
C'  is  the  line  parallel  to  £-axis  at  mid  radius.  The  calculation  points 
are  located  midway  between  the  singular  points  The  induced  velocities 
=re  calculated  at  points  very  near  the  blade  at  A0  = 2°  and  1/2,  1/4, 

1 8,  1/16-spacing  points  (See  Figure  50). 

It  should  be  pointed  out  that  in  the  present  model  every  blade  has 
the  same  distribution  of  and  without  regard  to  the  effects  of 
different  stagger  angles  (A)  and  the  number  of  the  blades  (NB) , which  is 
directly  related  to  solidity.  To  avoid  a misleading  conclusion,  the 
velocity  components  due  to  vortex  and  source  will  be  shown  separately  to 
pcxnt  out  the  effects  of  NB,  X and  V 
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The  typical  tangential  velocity  distribution  at  the  hub  induced  by 

NB -radial  vortex  lines  is  shown  in  Figures  51  and  52  for  the  same  hub- 

to-tip  ratio  v = 0,3  and  number  of  the  blades  NB  = 1 but  at  different 

stagger  angles  (A  = 0°  and  A = 54°)  to  show  the  effect  of  stagger  angle 

on  the  three-dimensional  solution  The  abscissa  represents  tangential 

direction  and  s is  equal  to  one  spacing  angle  (s  = 360°/NB) . Figure  53 

shows  the  ratio  u Q /u  » in  the  chordwise  direction  for  A = 0°, 

2d  3d 

A0  = 2U  at  the  hub,  mid  radius  and  tip.  The  angular  location  A0  = 2° 


is  chosen  since  the  maximum  discrepancy  in  the  tangential  velocity 

occurs  near  this  location.  As  mentioned  in  Section  4,3,1,  the  effects 

of  the  number  of  the  blades  on  the  tangential  velocity  are  quite 

similar  to  those  of  the  lifting  line  model.  Namely,  for  large  NB,  the 

ratio  u - /u  n goes  to  unity.  The  maximum  difference  between  the 
v62d  ve3d 

two  solutions  (about  30%)  occurs  when  NB  = 1,  v = 0,3,  The  ratio 

uv0  /uv0  dn  vicinity  a fading  edge  is  sometime  very  large, 

2d  3d 

but  the  absolute  value  of  u^  at  this  point  is  usually  small  due  to 

3d 

cancellation  of  the  induced  flow  due  to  F , and  T Hence,  the 

i,l  1,2 

ratio  u „ /u  p at  these  points  should  not  be  emphasized.  In  Figures 
V 2d  3d 

5a  and  55,  the  tangential  velocities  at  two  radial  locations  n = V (hub) 
and  n = 10  (tip)  for  NB  = 3,  A = 18°  and  v = 0-3  are  shown.  The 
difference  between  two-dimensional  and  three-dimensional  solutions  are 
larger  at  the  hub  than  at  the  tip  The  absolute  values  are  generally 
larger  at  the  hub  Generally,  the  difference  between  the  two  solutions 
a*  the  mid  radius  are  negligibly  small-  The  two-dimensional  values  are 
more  uniform  across  the  passage  than  the  three-dimensional  values  at  the 
hub  and  this  trend  is  reversed  at  the  tip  Two-dimensional  solutions 
underestimate  the  maximum  tangential  velocity  u^p  at  the  hub  and 


at  the  hub  and 


overestimated  the  value  at  the  tip.  The  effects  of  the  hub-to-tip 

ratio  v and  the  number  of  the  blades  (NB)  on  the  difference  between  two 

solutions  are  more  dominant  than  those  due  to  stagger  angle  X, 

The  effects  of  the  number  of  the  blades  on  the  induced  axial 

velocity  u , due  to  the  NB-senes  of  radial  vortex  lines  show  the  same 
v? 

trend  as  the  behavior  of  the  axial  velocity  in  the  "Lifting  Line  Model" 
mentioned  in  section  4.3  1,  Namely,  for  large  NB,  the  axial  velocity 
goes  to  zero  quite  rapidly,  especially  at  the  hub.  Therefore,  the  case 


of  NB=1  is  selected  to  discuss  the  effects  of  stagger  angle  X,  In 

Figures  56  and  57,  the  behavior  of  the  axial  velocity  v.  _ and  u 

Vo3d  ^2d 

are  shown  for  widely  varying  stagger  angle  X = 0 and  54  for  NP=1, 


V = 0.3  at  the  hub.  The  trend  is  almost  similar  for  both  cases,  but  the 


discrepancy  at  large  stagger  angles  is  slightly  larger  than  that  at 
small  stagger  angles.  To  show  these  effects  clearly,  the  axial  velocity 
ratio  u , /u  „ for  NB=1  at  L/ c cos  X = 0.0  and  n = v (hub)  are  shown 

V?2d  VSd 

m Figure  58.  The  axial  location  C,/ c cos  X = 0,0  is  selected  since 


maximum  difference  between  two  solutions  occurs  near  this  location.  The 


axial  velocity  ratio  is  larger  for  larger  stagger  angles  at  both  hub  and 
tip,  and  is  also  amplified  by  small  hub-to-tip  ratio.  Incidentally,  the 
ratio  at  large  axial  locations  should  not  be  emphasized  due  to  small 
absolute  values.  The  departure  at  the  hub  is  larger  than  those  at  the 
tip  and  at  mid  radius 

The  radial  velocities  induced  by  the  NB-senes  of  radial  vortex 
lines  are  discussed  below.  The  radial  velocity  distributions  u /vCTT” 

vn  ir  ^ 

is  plotted  along  the  tangential  direction  m Figure  59  at  the  mid  radius 
and  at  C/c  cos  ) - 00  for  X = 18°,  v = 0.6  and  0.3  for  various  values 


of  NB 


t is  clear  that  the  small  hub-to-tip  ratio  causes  larger  radial 


50 


velocity  thar  the  large  hub-to-tip  ratio,  and  the  larger  number  of  the 

blades  makes  the  radial  velocity  distribution  uniform  along  the 

tangential  direction  The  effect  of  the  stagger  angle  A is  relatively 

smaller  than  that  of  hub-to-tip  ratio  and  number  of  blades  (NB) . 

Generally,  the  order  of  magnitude  of  the  radial  velocity  induced  by 

NB-senes  of  radial  vortex  lines  is  approximately  15%  of  the  maximum 

tangential  and  axial  velocity  components  for  V = 0.6  and  30%  for  v = 03 

Unlike  the  vertex  system,  the  tangential  velocity  due  to  the  NB 

series  of  radial  source  lines  u„  decreases  relatively  slowly  along 

Sy3d 

the  axial  direction.  In  order  to  examine  the  effects  of  the  stagger 

angle  A,  two  chordwise  planes,  such  as  £/ c cos  A = 0 and  1.0,  are 

selected  The  ratio  u a /u  „ is  plotted  along  0 axis  for  v = 0.3  and 

s02d  S03d 

NB  - 1 at  the  hub  in  Figures  60  and  61.  Also,  the  effects  of  hub-tc- 

tip  ratio  (v)  cn  the  ratio  u „ /u  . for  v = 0.6,  NB  = 1 at  the  hub 

SH2d  S 3d 

are  plotted  in  the  same  figure  only  for  A = 0 and  54  The  effect  of 

the  hub-to-tip  ratio  v is  more  dominant  than  that  of  A. 

The  axial  velocity  due  to  the  NB -series  of  radial  source  lines 

u.  and  u is  usually  larger  than  the  corresponding  tangential 
'2d  5 "3d 

velocity  u „ and  u A , and  different  from  the  induced  velocity  due 
S02d  £03d 

to  a vortex  system  The  axial  velocity  u due  to  the  NB-series  of 

5 -3d 

radial  source  and  sink  lines  causes  quite  a large  discrepancy  from  the 

two-dimers lonal  solution  for  large  NB,  especially  at  the  hub  for  small 

hub-to-tip  ratio  as  mentioned  in  the  "Source  Line  Model".  The  effects 

of  the  NB  on  the  ratio  (u  /u  ) for  zero  stagger  angle  (A  = 0)  are 

S;2d  S'3d 

shown  in  Figure  62,  (NB  = 1,  3,  5,  20  for  V = 0 3 and  N'B  = 1,  20  for  v = 
06’  For  reference,  the  limiting  values  of  (u  /u  „ ) at  far  down- 

«2d  sSd 

stream  £ = <*  are  also  shown  in  the  graph  They  are  (lrv/2ri)  = 2,17  at 


r 


^ ■ 
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rhe  hub  and  0,65  at  the  tip  for  v = 0,6,  At  the  point  close  to  the 

leading  edge  and  trailing  edge,  the  ratio  has  no  meaning,  since  the 

axial  velocity  u is  very  small  The  ratio  (u  /u  ) at  the  mid 
-3d  b 2d  s 3d 

radius  is  very  close  to  unity  except  near  the  leading  edge  and  the 

trailing  edge  The  effects  of  stagger  angle  on  the  ratio  (u  /u  ) 

s;2d 4  S '3d 

at  ’"he  hub  for  NB  = 20  are  shown  in  Figure  63  in  the  region  very  close 
to  the  singular  point  (A0  =•  2°)  . The  effects  of  stagger  are  negligible 
for  large  hub-to-tip  ratio  (v  = 0,6),  and  the  effects  of  NB  and  v are 
more  dominant  than  those  of  stagger. 


The  radial  velocity  distribution  u //_.  2 induced  by  NB  series  of 

sr  IQ 

J 

radial  source  lines  are  shown  in  Figure  64  along  the  tangential  direction 
for  A = 18°  at  £/c  cos  A - 0 at  mid  radius  for  several  NB  and  V, 

Similar  to  the  vortex  system,  the  radial  velocities  due  to  the  NB  series 
of  radial  source  lines  for  small  hub-to-tip  ratio  (v  = 0,3)  are  larger 
than  those  of  large  hub-to-tip  ratio,  The  larger  number  of  the  blades 
NB  makes  ’-he  distribution  of  the  radial  velocity  uniform  along  the 
tangential  direction  However,  the  radial  velocity  induced  by  the  NB 
series  cf  radial  source  lines  for  large  NB  is  larger  than  that  for  small 
NB  Therefore,  for  the  case  of  large  NB  and  small  hub-to-tip  ratio,  the 
radial  velocity  due  *o  source  lines  is  not  necessarily  small  compared 
with  other  velocity  components.  The  effect  of  the  stagger  angle  A is 
relatively  smaller  than  those  due  to  NB  and  v. 

4 ? . 4 Effects  of  the  Radial  Source  Lines 
of  Variable  Strength 

All  Tt>e  example  calculations  given  in  the  preceding  section  were 
earned  cuf  for  a series  of  radial  vortex  and  source  lines  of  constant 


snength  However,  it  is  not  necessary  for  a radial  source  line  to  be 


52 


of  constant  strength  as  mentioned  in  Section  2.2,  It  is  expected  that 
a series  of  source  lines  of  variable  strength  can  have  three-dimensional 
effects  directly.  In  order  to  examine  the  effects  of  the  radial  source 
line  of  variable  strength  along  the  radial  direction,  a three-bladed 
system  of  symmetric  profiles  ana  of  zero  stagger  angle  was  considered 
for  hub-to-tip  ratio  v = 0,6.  Each  blade  is  replaced  only  by  a pair  of 
radial  source  lines  Q.  and  Q,  The  locations  of  Q.  1 and  Q,  are 
shown  in  Figure  65,  and  the  values  of  them  are  given  by 

l = - 0,5n  + 0,9  (89) 

and 

Q = 0.5n  - 0.9,  (90) 

i > *• 

where  i = 1,  2 and  3,  For  the  sake  of  convenience,  the  inlet  flow 
velocity  on  which  the  three-dimensional  perturbed  velocities  are  super- 
imposed is  taken  to  be  unity.  For  several  angular  locations  the 
numerical  results  of  the  pressure  distribution  based  on  the  inlet  flow 
velocity  are  shown  m Figures  66  and  67  at  5 = 0 and  L = 0,5, 
respectively.  The  point  £ = 0 is  near  the  leading  edge  and  £ = 1.0  is 
near  the  trailing  edge  The  flow  field  corresponding  to  £ = 10  is  the 
same  as  rhat  at  5 = 0,  and  the  values  for  negative  9 are  the  same  as 
rhose  for  positive  9 by  symmetry.  At  £ = 0,  the  discrepancies  between 
the  two-  and  three-dimensional  analyses  are  quite  small  for  any  9,  and 
the  values  based  on  the  analyses  are  coincident  at  mid  radius  p = 08 
(see  Figure  66),  However,  in  the  region  near  the  suction  surface 
of  the  blade  at  ? = 0.5,  the  pressure  coefficients  based  on  the  three- 
dimensional  analysis  are  different  near  the  hub  and  tip  (see  Figure  67) . 
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system  with  constant  strength  of  the  radial  source  lines,  j = 0,2 

and  Q = - 0.2,  is  shown  m Figures  68  and  69,  The  strength  of  the 
i * ^ 

source  lines  is  taken  to  be  the  integrated  mean  value  of  the  preceding 
model-  Comparing  the  two  cases,  the  discrepancies  between  the  two- 
dimensional  unrolled  solutions  and  the  three-dimensional  solutions  are 
relatively  smaller  than  those  of  variable  strength  source  lines  Thus, 
the  effect  of  a source  line  of  variable  strength  upon  the  three- 
dimensional  flow  field  is  large  and  the  discrepancies  between  the  two- 
and  three-dimensional  analyses  are  amplified  by  the  variation  of  the 
strength  of  the  source  line  along  the  radial  direction  Furthermore, 
the  differences  between  the  pressure  coefficients  of  the  two  analyses 
vanish  at  mid  radius  for  the  hub-to-tip  ratio  v = 0 6 for  both  constant 
and  variable  strength  source  lines. 
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CHAPTER  V 
CONCLUSIONS 


The  flow  field  in  the  axial  flow  turbomachinery  encountered  in 
practice  is  three-dimensional  and  is  affected  by  viscosity  and 
compressibility.  Therefore,  the  three-dimensional  potential  flow 
theory  can  be  used  only  for  restricted  three-dimensional  flow  problems 
because  of  the  assumed  incompressible,  irrotational , and  inviscid  flow 
conditions.  The  three-dimensional  potential  flow  theory  can  play  an 
important  part  in  the  clarification  of  the  three-dimensional  effects  on 
the  flow  field  of  an  incompressible  perfect  fluid  The  three- 
dimensional  potential  flow  theory  based  on  the  analysis  of  Rossow  (3) 
and  Tyson  (2)  was  modified  in  order  to  predict  the  three-dimensional 
flow  effects.  The  blade  row  is  replaced  by  a series  of  radial  vortex 
lines  of  constant  strength  along  the  ladial  direction  and  a series  of 
radial  source  lines  of  variable  strength  along  the  radial  direction. 

The  perturbed  velocity  field  due  to  a single  radial  vortex  and  source 
line  of  constant-  strength  was  calculated  for  hub-ro-tip  ratios  v = 0,6 
and  0.3  The  perturbed  velocity  field  due  to  a radial  vortex  line  of 
constant  strength  for  v = 06  was  coincident  with  the  results  of 
Reference  (4),  For  the  perturbed  flow  field  due  to  a radial  source 
line,  there  is  no  other  numerical  result  available  for  comparison  The 
qualitative  behavior,  however,  is  very  similar  to  that  of  the  two- 
dimensional  unrolled  solutions  except  for  the  radial  velocity  components 
The  characteristic  behavior  of  the  three-dimensional  perturbed  flow 
field  due  to  a single  radial  vortex  and  source  line  of  strength  unity 
discussed  in  Sections  422  and  4.2,3  are  summarized  below 
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a)  The  qualitative  behavior  of  the  tangential  and  axial  velocity 
components  is  similar  to  that  for  two-dimensioral  unrolled  solutions. 

The  maximum  discrepancies  between  both  analyses  are  25%  for  v = 0 6 and 
30%  for  v = 0-3. 

b)  The  tangential  and  axial  velocity  components  due  to  a single 
radial  vortex  line  of  constant  strength  are  very  similar  to  the  axial 
and  tangential  velocity  components  due  to  a single  radial  source  line  of 
strength  unity  as  in  the  two-dimensional  unrolled  solution.  Also,  the 
behavior  of  the  velocity  components  due  to  a source  line  relative  to 
the  two-dimensional  values  is  always  opposite  to  that  due  a vortex  line. 

c)  The  behavior  of  the  radial  velocity  components  due  to  a single 
radial  vortex  and  a source  line  of  constant  strength  are  very  different 

d)  At  mid  radius  (q  = l+v/2) , the  tangential  and  axial  velocity 
components  of  the  three-dimensional  flow  field  are  very  close  to  those 
for  two-dimensional  theory  for  hub -to -tip  ratio  v = 06,  but  small 
discrepancies  between  the  two  analyses  were  found  for  v = 03 

e)  The  three-dimensional  effects  are  amplified  by  a small  hub-to- 
t.ip  ratio 

f)  The  discrepancy  between  two-dimensional  and  three-dimensional 

I 

solution  at  the  hub  (q  = v)  are  opposite  to  those  at  the  tip  (q  = 1 0) 

and  these  are  amplified  by  a small  hub-to-tip  ratio  for  both  a vortex 

line  and  a source  line.  • 1 

The  three-dimensional  potential  flow  theory  was  modified  to  solve 


the  three-dimensional  flow  field  around  a blade  row  of  NB  blades,  where 
each  blade  is  replaced  by  a series  of  radial  vortex  lines  of  constant 
strength  and  a series  of  radial  source  lines  of  variable  strength 
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Example  calculations  were  carried  out  for  the  "Lifting  Line  Model", 
"Source  Line  Model",  "Lifting  Surface  Model  Including  Thickness  Effects" 
and  "Effects  of  the  Radial  Source  Line  of  Variable  Strength",  In  the 
first  and  second  example  calcux  ions,  each  blade  is  represented  by  a 
single  radial  vortex  line  (T  = - TO)  or  a source  line  (Q  = 1,0),  and 
the  effects  of  the  number  of  the  blades  NB  and  hub-to-tip  ratios  are 
examined,  The  effects  of  the  stagger  angle  A together  with  NB  and  v 
are  examined  in  the  third  example  calculation  in  which  each  blade  is 
replaced  by  a series  of  radial  vortex  lines  and  source  lines  of 
constant  streng+h  , The  calculations  are  carried  out  for  stagger  angles 
0°,  18°,  36°  54°,  number  of  the  blades  NB  = 1,  3,  5,  20  and  hub-to-tip 

ratio  v = 0 6 and  0,3-  Finally,  the  effects  of  the  radial  source  line 
of  variable  strength  are  investigated  fox  simple  models  The 
characteristic  features  of  the  three-dimensional  flow  field  are 
summarized  below. 

a)  The  effects  of  the  small  hub -to -tip  ratio  on  the  discrepancy 
between  the  two-  and  three-dimensional  solutions  are  greater  than  those 
of  the  number  of  the  blades  NB  and  the  stagger  angle  A. 

b)  Effects  of  the  stagger  angle  A on  the  discrepancies  between  two 
solutions  is  usually  smaller  than  those  due  to  NB  and  v, 

c)  Small  number  of  blades  (NB)  cause  larger  discrepancies  between 
two  solutions,  especially  at  the  hub,  except  for  the  axial  and  radial 
velocity  components  due  to  NB  series  of  radial  source  lines.  Alterna- 
tively, the  axial  and  radial  velocity  components  due  to  source  lines 
are  amplified  by  large  number  of  the  blades  (NB)  together  with  small  v. 
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d)  The  radial  source  line,  of  variable  strength  causes  large 
discrepancy  between  two  solutions  even  for  relatively  large  hub-to-tip 
ratio,  but  at  the  mid  radius,  the  discrepancies  vanish  for  V = 0.6 
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APPENDIX  A 


Specif ica^i on  of  the  Boundary  Condition  for 
Radial  Vortex  Line  and  Source  Line 


Specif icat ion  of  the  Boundary  Conditions  for  a Radial  Vortex  Line 

The  single  radial  vortex  line  of  constant  strength  I located  at 
c = 0,  0 = 0°  is  split  into  ydto  (see  Figure  70)  and  distributed  over 
- 6/2  < 8 < 6/2  as  given  by 


6/2 

/ ydoo  = y'6  = T - 


The  induced  velocity  Aq  at  (p,  6,  6/2)  due  to  ydto  is  given  by 


where 


Aq  - ■yda'/Arh  (cos  a - cos  , 


n cos  (0  - to)  - m 


, (A3) 


p + m4  * (Ac/2)  - 2mn  cos  (0  - to) 


AL  _ l t m - p cos  (8  - to) 

P L p4  * (Ji-^m)  ‘ + (AC/ 2) -2  (Jt-m)pcos  (0  -to) 


sin  a = 


Ac  /2 


p1  sin4  (8  - to)  + (A£/2) 


Therefore,  the  tangential  velocity  component  of  the  induced  velocity  due 
to  the  total  of  the  split  radial  vortex  lines  is  given  by 
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Integrate  and  take  a limit  as  A£  * 0;  then, 

Mm  vc  = p(n)/2n,  for  - 6/2  < 9 < 6/2  (All) 

A£-K) 
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APPENDIX  B 

Elimination  of  the  Summation  About  1 in  the  Series  Summation 


Usually  the  blade  row  in  turbomachmeiy  is  constructed  by 

equally  spaced  blades  of  the  same  profile.  Therefore,  the  distribution 

of  the  series  of  vortex  and  source  lines  is  the  same  for  any  blade. 

Now,  0 is  the  angular  coordinate  of  the  location  of  the  l-th  vortex 

i,l  J 

or  source  line  for  the  i-th  blade.  From  the  periodic  characteristics. 


we  get 


0 =6,  + ( 2 n/NB ")  (i  - 1)  , 

i»J  1,1 


Using  Equation  (A12),  we  can  eliminate  the  summation  with  respect  to  i 
in  Equation  (39)  in  Section  3,1,  Thus, 


sin  m(0-0  .)=  sin.  m(0-0  .)cos  2miT(i -1) /NB-cos  m(0-0.  ,)sin  2miT(i-l)/NB 

1 > 1 1,1  1 , 1 

(A13) 


>(l,m)  = I exp  (-K  (£,-£.))  • F (K  V)  • Z (K  n)  (A14) 

’ - r + mn  j v ^ mn,  nr  mn  J 


JVS  r NB 

0 = - Z r-  l (0-0, 

v * , 4ir  , li 

j=l  i=l 


•If11-1’*2  2 "(9-VC0ST^i> 

m=  1 J 


- cos  m (0-0j j ) sin  — ] S(j,m)) 


■ 
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JVS 

y 


NB 


l ii(NB(0-9  )-(NB-1)tt+2  Uysm  m(0-0.  H C0S£iEiilLL 

j-1  13  m- 1 a 53  NB 


1=] 


NB  . . 

cos  m (0-0lj)  Z sin  ] S(j,m)) 


(A15) 


and 


NB  2mT>(i-n 

I cos  — n-  -g-- = cos  ((1-1/NB)  • imr)  »sin  mn/sin (mn/NB) , 

i = l 

= 0 

where  m f l • NB , When  m is  a multiple  of  NB  Cm  = £ • NB) , we  get 


NB  NB 

I cos  — m-gg1 — - =E  cos  (i-1)  *2  Air)  = NB 


i=l 


i = l 


Simi  !s.r  ly , 


NB 

I sin  = sin  ( [1  -l/NB)rMT)  *sin  m rr/sin (mr/NB) 

1 = 1 


= 0 


where  m t £ • NB,  and  when  m ~ l • NB, 


NB  2mrfi  _n  NB 

E sin  — - E sm  (i  - 1)  * 2JZtt)  = 0 
i=l  1 r 1 


Therefore,  the  velocity  pofential  reduces  to 


JVS 


<t>  = 2 NB  E (6-0.  -(1  - 1/NB)tt+2  E 


00  sin  k (0-0.  .) 


3 = 1 * ' 13 £,=  1 k 


ill  . 


S(j.k)), 


(A16) 

where  k = £ • NB  = 1 , 2 , 3 , . . . In  the  same  way,  the  summation  with 


respect  to  i in  equation  (46)  can  be  eliminated 
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APPENDIX  C 

Related  Bessel  Function  and  Characteristics  of  the  Series  Expression 


a)  Function  Z (K  i'll 
m mn 


Z (k  rj  = J (K  n) 
m • mn  • m mn 


3Y  (K  n) 
m mn 

3n 


- I 


n=i 


9J  ( k n) 

m , — — ] y (K  n) 

id  , m mn 

n=i 


(A17) 


From  Equation  (A17) 


z ;k  n)  = k t J (x) • t 

m mn  mn  m A 3X 


X=K 


9J  (X) 

i-fr-1  • Y» <x))  - 


mn 


X=K 


mn 


where  y = K r Then  by  Lommel's  formula  for  Bessel  Function, 
A mn 


[Z  (K  n)l  = K [— ] - - 

L m mn  u 1 mn  ry  * 


(A18  i 


and  also  frcm  ehe  bcundary  condition  given  by  Equation  (11)  in  Section 
2 1, 


J ' (K  u)  J ’ . K ) 

m mn  m mn 


( A 1 9 4 


Then , 


Y '(K  v>)  r ' ( K ) 

m mn  m mn 


Y * (K  ) 

Z iK  n)  - -2 52—  Z K n)  . 

m mn  x,  , ,,,  . , m mn 


(A20j 


where 


Y '(K  v) 
m v mn 


Z K - J (K  n)  Y ' (K  v)  - J ’ (K  v)  Y (K  n) 
m mr  ' m mn  m mn  m mn  m mn 


Therefore 


ayx)  3Jm(X) 


(Z  (K  rlj  = (Z  (x)J  * K [J  (x)  3x  3x  «n 
m mn  m mn  m 


n v 


X=k  v 
A mn 


X=vK 


mn 
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by  Lommel's  formula 


[Z  (K  r)l  - K [2/tty]  = 2/vtt 

mv  mn  ' mn. 

X — In  V 

n=v  A mn 


Hence , 


0 Y '(K  ) 

fZ  (K  n)  1 = — 
m mn  vtt 


Y ' (K  v) 
n=v  m mn 


(A21) 


bl  Integral  and  Derivatives  of  Z (K  p) 
6 m mn 

From  boundary  condition  (A20) , 


9Z  (K  n) 
m mn 


dn 


= 0 


rpl  ,v 


(A22) 


3 Z (K  rv 

m irn 

N 2 

on. 


2 2 2 

— (m  -K  ) when  n = 1 
TT  v mn 


— (^=-  - K 2 
vtt  ^2  mn. 


Y ’ (K  ) 
m mn 

Y ’ (K  v) 
m mn 


when  n = v 


( A23 1 


and 

1 


_ 2 „ , . 1 mp__2_.,„ 

r,z  ,K  Pi  dp  = — -=■  SZ  (SldS  = 

m mn  ^ 2 m 

v mn  K v 

mn 


((1  - siz  2(S)*Z  2(S))) 

,,  2 m m K v 
S mn 


(A24) 


where  I Z ' i S) ] 

m S- K , K v 
mn  mn 


= 0 


c)  Closed  Form  for  Series  Expression 


By  Dim's  expansion, 
1 


<*  rZ  (K  p)d.P 

v m mn  _ . , 

I — ; Z tK  P)  = 1, 

.1  - nr  mn 

/ nZ  (K  p)dn 
1 m v mn 
v 


CA25) 
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using  Equation  (A18) , we  obtain 


1 

00  n)dn 

^ v m mn 


n=l 


1 2 

j riZ  (K  r|)dri 
m mn 


7T 

? ° 


(A26) 


The  tangential  and  radial  velocity  components  due  to  a single  radial 
vortex  line  and  also  the  axial  velocity  component  due  to  a single  radial 
source  line  are  zero  at  £ = 0, 


{uv9} 


1 

r 00  “ / nz  (K  n)dn 

: (1  + 2 E cos  m e E ~ — z (K  n))  = o , 

4r|T  v , .1  - m mn  * 

£=0  ! n Z (K  p)dn  (A27) 

J m mn 


(u 


sm  m 9 


2r 


CO 


m=l 


n=l 


( r\Z  (K  n)dr| 
'V  m mn 

1 2 

/ r\Z  (K  n)dr. 
‘ m mn 

v 


8Z  (K  n) 
m mn 

3n 


= 0 (A28) 


and 


1 2 

tUS£;  * 4-  2 

4=0  1_V 


1 

1 - L QC0)Z  (K  n)dn 

QCn)dn±  £ 7^ VKon^ 

V n J nZo2(Konn)dn 


V 


00  00  / q (n)  z (K  p)dn 

± 2 E cos  m 9 I -y 2L-J52 Z^K^n))  = ° (A29) 

mil  n=1  ( nz  2(k  n)dn 

m v mn 

v 

The  tangential  and  radial  velocity  components  due  to  a single  radial 
vortex  line  of  constant  strength  should  have  an  inflection  point  at 
4-0,  then, 
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~2 
a U. 


V0, 


. _r_ 

'+  2-rrn 


£ cos  m 0 £ K 


2 fv 


t=0 


m=l 


n=l 


mn 


/ nz  2(k  n)dn 

; v m mn 


z (K  n)  = o 

mv  mn  J 


(A30) 


and 


32u 

[— 


H 


5=0 


« 00  / nzm(K  n)dn 

. - r , sinme  r „ 2 v m mn 

' « » Zj  ~ Zj  k — — 

2tt  , m _ mn  1 ^ 3p 


9Z  (K  n) 

nK  mn  ’ 


m=l 


n=l 


/ nz  (K  n)dn 

1 m v mn  J 


= 0 . 
(A31) 


Also,  the  axial  velocity  component  due  to  a single  radial  vortex  line 
of  constant  strength  has  a maximum  value  at  £ = 0;  then. 


3u  r 00  • a “ J nZm(Kmnn)dn 

[ -Jy  ■ ' ] = zJ-  £ — n-m-9  £ K 2 -Z 

5=0  m 1 n 1 / nz  2(K  n)dn 

1 m mn  J 


•Z  (K  ) = 0 
m mn 


v 


(A32) 


Similarly,  for  the  three-dimensional  flow  field  due  to  a single  radial 
source  line  of  variable  strength,  the  tangential  and  radial  velocity 
components  have  a maximum  value  at  £ = 0;  then. 


/1Q(n) Z (K  n) dp 


3uCQ  . 00  «>  ' ■'  mv  mn 

f Ot/  -j  « 1 n • r\  T' 

] = * JnrT  E1m  Sin  m 0 E,  — : 

5=°  m 1 n_1  /nz  (K  Il)dri 

' m mn 


•Z  (K  n)dn  = 0 (A33) 

nr  mn  J K 


V 


and 
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3u  , . / Q(n)Zo(Konn)dn  3Z 

. sn , _ 1 r-  y 0 on 

■3?  4tt  ■ 7 1 2 3rl 

• c«o  n-'  / nz  "(K  n)dn 


I Q(n)Z  (K  n)dn  JZ  (K  n) 
oo  °°'  m mn  mmn 

i 2 I cos  m 6 £ -j ) = 0 . 

m'1  n=1  / n.z  2tk  n)dn 

J m *■  mn 
v 

The  axial  velocity  component  has  an  inflection  point  at  £ = 0;  then, 


2 < Q(n)Z  (K  n)dn 

) u , « 1 o on 

— —}  = — r-  e k -y z (K  n) 

2 4tt  *•  f 9 on  1 7 o on 

?c  t=o  n f nz  (K  n)dn 

- o on 

v 


/ Q(n)Z  (K  r)dn 
oo  oc  j w u m mn 

2 I cos  m 6 l K 2 (K  n))  = 0 

, , mn  1 _ m mn 

m'  1 n=1  / nZ  2(K  n)dn 

* m mn 


69 


APPENDIX  D 

Approximate  Formula  for  the  Solution  of  the 
Three-Dimensional  Flow  Field 


The  solutions  of  the  three-dimensional  flow  field  due  to  a single 
radial  vortex  line  and  source  line  include  series  summations  of  Bessel 
functions  and  trigonometric  functions  and  are  quite  complicated,  A 
modification  of  the  three-dimensional  . olutions  will  be  attempted  here 
The  qualitative  behavior  of  the  three-dimensional  flow  field  is  similar 
to  the  two-dimensional  unrolled  solution,  The  perturbed  velocities  are 
approximately  identical  for  both  two-dimensional  and  three-dimensional 
cases  at  mid  radius  for  any  hub-to-tip  ratio.  Rossow  (3)  suggested  the 
use  of  Taylor  expansion  for  the  exact  solution  of  the  three- 
dimensional  flow  field  due  to  a single  radial  vortex  line  as  an  approxi- 
mate formula  valid  for  hub-to-tip  ratios  close  to  unity.  However,  this 
approximate  formula  may  not  be  valid  for  small  hub-to-tip  ratios 
because  of  the  assumption  made  in  the  process  of  the  Taylor  expansion 
In  the  process  of  calculating  the  three-dimensional  flow  field,  dupli- 
cate calculations  are  required  in  order  to  superimpose  the  effects  of 
the  vortex  and  source  lines.  Once  the  basic  solutions  are  obtained,  the 
superposition  of  the  flow  field  for  several  vortex  and  source  lines 
involves  only  simple  algebraic  equations.  Therefore,  if  the  basic 
solutions  are  represented  by  simple  approximate  formula  based  on 
numerical  results  of  the  basic  solutions,  these  approximate  formula 
are  valid  for  any  calculation  for  the  same  hub-to-tip  ratio.  Several 
approximate  formulations  are  suggested  below: 


a)  Approximate  velocity  potential  and  velocity  components  due  to 
a single  radial  vortex  line  of  constant  strength 


. r . -1  tan  0/2 

($J  = “ 2?  an  

ap  tanh  B£ 


(A36) 


F sinh  2 

4irr|  cosh  2BC  - cos  0 ’ 


. _ _T_  B'  C sin  6 

Uvri ^ 2tt  cosh  2B£  - cos  0 
ap 


. _ _T_  B sin  0 

■vt,  2tt  cosh  2B£  - cos  0 5 
ap 


where  B - B(rj) . B'  = 3B(p)/9r),  These  are  similar  to  the  two- 
dimensional  unrolled  solutions  except  for  the  radial  velocity  (see 
Appendix  E) . The  functional  behavior  of  the  approximate  formula  are 
found  to  be  quite  similar  to  the  three-dimensional  exact  solutions.  The 
approximate  tangential  and  radial  velocity  components  have  the  same 
limiting  values  at  C -*■  ± 0 and  C -*■  ± 00  as  the  three-dimensional  exact 
solutions  Also,  the  approximate  tangential  and  radial  velocity 
components  have  inflection  points  at  £ = 0 and  the  axial  velocity 
component  has  a maximum  at  £ = 0,  However,  the  functions  B(t|)  and  B'(rO 
are  unknown  Bur,  from  numerical  results  of  the  exact  solution  at  the 
particular  r and  0,  the  values  B(n)  and  B'(n)  in  Equations  (A37) , (A38) 
a™ d (A39)  can  be  determined  The  development  of  these  approximate 
solutions  ?.s  planned  for  the  future. 
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b)  Approximate  velocity  potential  and  velocity  components  due  to 
a single  radial  source  line  of  constant  strength: 

The  approximate  formula  for  the  three-dimensional  flow  field  due 
to  a single  radial  source  line  of  variable  strength  can  be  written  as 


C*S>  * 4,(1?V)C  ^ <2<C°Sh  2C?  - C°S  e>>  ’ 

ap 


(A40) 


sin  6 


S0  4TT(l+v)Cn  cosh  2C^  - cos  0 ’ 

aP  ’ 


(A41) 


(uc  ) * - ■ (1  + £n(  2 (cosh  20;  - cos  0))  - 2 2Cg 

Srr  47T(lrv)C  v C ^ cosh  2CC  - cos0 

ap  (A42) 


and 


(ucr)  = 


sin  2Ct; 


S^ap  2tt(1+v)  cosh  2C£  - cos0  ’ 


(A43) 


where  C = C (n) , C'  = 9C(r|)/9n,  and  C = l/2n  correspond  to  the  two- 
dimensional  unrolled  solutions  (see  Appendix  E) , Again,  the  functions 
C(r,1  and  C'(n)  are  unknown,  similar  to  the  vortex  system.  These 
quantities  will  be  determined  for  fixed  q and  0.  The  values  C(q)  and 
C'fq)  so  derived  are  valid  for  the  estimation  of  the  velocity 
components  along  the  t;  axis  for  particular  n,  0,  v.  Values  of  B,  B', 
and  C'  of  the  present  formula,  whose  functional  behavior  are  quite 
similar  to  the  corresponding  exact  solutions,  will  be  determined  by 
using  the  exact  numerical  results  at  a few  points  for  particular  n,  0 
and  t,  Then,  the  approximate  velocity  components  at  any  £ coordinate 
can  be  easily  determined 
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APPENDIX  E 


Two-Dimensional  Unrolled  Solution 


The  velocity  potential  due  to  a vortex  of  strength  T and  a source 
of  strength  Q located  at  £ = 0 and  0=0  are  given  by. 


. T . -1  tan  0/2 

$ = - x—  tan  - — ■■■■  — 

v 2tt  tanh  ^/2r\ 


(A44) 


$s  = S-  iln  (2(cosh  £/q  - cos0)) 


The  corresponding  velocity  components  are: 

T sinh  £/r| 

Uv0  47TT1  cosh  C/n  - cos  0’ 


(A46) 


_T sin  0 

JvC  - 4Tm  cosh  £/r) 


cos  0 ’ 


Q sin  0 

1SQ  4iTri  cosh  ?/n  - COS0 


(A48) 


Q sinh  C/n 

4Ttn  cosh  ^/n  - cos0  ’ 


(A49) 


where  the  positive  direction  of  F is  based  on  the  right-hand  screw  rule, 
and  the  positive  and  negative  values  of  Q correspond  to  a source  and 
sink,  respectively.  The  tangential  and  axial  velocity  components  due 
to  a series  of  radial  vortex  lines  of  constant  strength  and  due  to  a 
series  of  radial  source  lines  of  variable  strength  replacing  the  equally 
spaced  blade  rows  of  NB  blades  are  given  by, 
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Figure  11:  Tangential  Velocity  Distribution  ur,  at  the  Hub( f = 0.6)  due  to  a Single  Radial 
Vortex  Line  of  Strength  Unity  for^=  0.6 


Figure  12:  Tangential  Velocity  Distribution  u^  at  the  Mid  Radius(?  - 0.8)  due  to  a Single  Radial 
Vortex  Line  of  Strength  Unity  for  p = 0.6 
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Axial  Velocity  Distribution  u yj  at  the  Tip(  = 1.0)  due  to  a Single  Radial  Vortex 
of  Strength  Unity  for  ^ = 0.6 


Figure  17:  Radial  Velocity  Distribution  u at  the  Mid  Radius(£  = 0.8)  due  to  a Single  Radial 
Vortex  Line  of  Strength  Unity  forp=  0.6 


Tangential  Velocity  Distribution  usa  at  the  Hub( 7 = 0.6)  due  to  a Single  Radial  Source 


Figure  19:  Tangential  Velocity  Distribution  u,^  at  the  Mid  Radius(,f=  0.8)  due  to  a Single  Radial 


of  Strength  Unity  for i)=  0.6 


tribution  at  the  Mid  Radi  us  (?  = 0.8)  due  tc  a Single  Radial  Sourc. 


Figure  2U:  Radial  Velocity  Distribution  us^  at  the  Mid  Radius(£*  0.8)  due  to  a Single  Radial 
Source  Line  of  Strength  Unity  for*>*  0.6 
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25:  Tangential  Velocity  Distribution  at  the  Hub(£*  0.3)  due  to  a Single  Radial  Vortex 
Line  of  Strength  Unity  forV=  0.3 


Tangential  Velocity  Distribution  uw  at  the  Mid  Radiu3(  1 = 0.65)due  to  a Single  Radial 
Vortex  Line  of  Strength  Unity  fory=0.3 


28:  Axial  Velocity  Distribution  ur^ at  the  Hub( 7 = 0.3)  due  to  a Single  Radial  Vortex  Line 
of  strength  Unity  for*^=  0.3 
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al  Velocity  Distribution  u,y  at  the  Tip( f = 1.0)  due  to  a Single  Radial  Vortex  Line 
Strength  Unity  for  \)  = 0.3 


Radial  Velocity  Distribution  uK^  at  the  Mid  Radlus(f*  0.65)  due  to  a Single  Radial 
Vortex  Line  of  Strength  Unity  for  V = 0.3 


Figure  32:  Tangential  Velocity  Distribution  us9  at  the  Hub(  f = 0.3)  due  to  a Single  Radial  Source 
Line  of  Strength  Unity  forv' = 0.3 
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Tangential  Velocity  Distribution  at  the  Mid  Radiu3(  J a 0.65)  due  to  a Single  Radial 
Source  Line  of  Strength  Unity  for  V - 0.3 


Tangential  Velocity  Distribution  us^  at  the  Tip(  'J  =*  1.0)  due  to  a Single  Radial  Source 
Line  of  Strength  Unity  for^  = 0.3 


al  Velocity  Distribution  u5j.  at  the  Hub(  1 = 0.3)  due  to  a Single  Radial  Source  Line 
Strength  Unity  forv  =0.3 
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Figure  42:  Ratio  of  the  Axial  Velocities  induced  by  the  Radial  Vortex 
Lines  of  JVS=1  at  5 »0.1,?  and 
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Radial  Velocities  induced  by  the  Radial  Source 


Fig.  1*9:  Singular  points  and  calculation  points  on  -plana 


Fig. $0:  Singular  points  and  calculation  points  on  (y ' ,g  )-plane 
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Distribution  of  the  Tangential  Velocities  at  the  Hub  Induced  by  HB  Series 
of  Radial  Vortex  Lines  of  JVS-8  for7\-0°.  NB-1  and  ^-0.3 


52:  Distribution  of  the  Tangential  Velocities  at  the  Hub  Induced  by  KB  Series 
of  Radial  Vortex  Lines  of  JVS“8  for7\»5U°,  KB-1  and  *>■(). 3 


Distribution  of  the  Tangential  Velocities  at  the  Tip  induced 
Radial  Vortex  Lines  of  JVS= 8 forA  = l8°,  NB=3  andP=0.3 


Fig. 56:  Distribution  of  the  Axial  Velocities  at  the  Hub  induced  by  NB  Series 
Radial  Vortex  Lines  of  JVS=8  forA=0°,  NB=1  andi^=0.3 


Distribution  of  the  Axial  Velocities  at  the  Hub  induced  by  NB  Series  of 
Radial  Vortex  Lines  of  JVS=8  for  a =514°,  NB=1  and  ^=0.3 
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Figure65  : Distribution  of  the  Source  of  Variable  Strength  and 
Coordinate  of  the  Singular  Points  ar.d  Calculation 
Points 
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Figu re66  : Pressure  Distribution  due  to  Source  Lines  of  Variable  Strength 
atJ^O  for  !I3=3  and  ^=0.6 
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Figure67  : Pressure  Distribution  due  to  Source  Lines  of  Variable  Strength 
at  J=0.5  for  N3=3  ar.d  i>~0.6 
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turbed flow  fields  are  solutions  of  the  linear  Laplace  equation  and  the  boundary 
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radial  variation  of  the  strength  of  the  radial  source  line  were  examined.  The  three 
dimensional  effects  are  found  to  be  appreciable  for  low  hub  tip  configuration  with 
fewer  number  of  blades. 
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